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ABSTRACT 

Using new spectra of 88 K and M-type subdwarfs, we consider novel methods 
for assigning their spectral types and take steps toward developing a comprehen- 
sive spectral sequence for subdwarf types K3.0 to M6.0. The types are assigned 
based on the overall morphology of spectra covering 6OOOA to 9OOOA. The types 
and sequence presented link the spectral types of cool subdwarfs to their main 
sequence counterparts, with emphasis on the relatively opacity-free region from 
8200-9000A. When available, supporting abundance, kinematic, and trigonomet- 
ric parallax information is used to provide more complete portraits of the observed 
subdwarfs. We find that the CaHn (n = 1-3) and Ti05 indices often used for 
subdwarf spectral typing are affected in complicated ways by combinations of 
subdwarfs' temperatures, metallicities, and gravities, and we use model grids to 
evaluate the trends in all three parameters. Because of the complex interplay of 
these three characteristics, it is not possible to identify a star as an "extreme" 
subdwarf simply based on very low metallicity, and we suggest that the modi- 
fiers "extreme" or "ultra" only outline locations on spectroscopic indices plots, 
and should not be used to imply low or very low metallicity stars. In addition, 
we propose that "VI" be used to identify a star as a subdwarf, rather than the 
confusing "sd" prefix, which is also used for hot O and B subdwarfs that are 
unrelated to the cool subdwarfs discussed in this paper. 

Subject headings: stars: subdwarfs — stars: abundances — stars: late-type — 
stars: classification — stars: fundamental parameters — stars: temperatures 
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Introduction 



The HR diagram is the most important map of stellar astronomy. It provides a relatively 
straightforward method for separating different stellar luminosity classes, e.g. supergiants, 
bright giants, giants, subgiants, main sequence dwarfs, and white dwarfs, using their colors 
and luminosities. However, spectroscopic and trigonometric parallax results have revealed a 
seventh distinct stellar luminosity class — the subdwarfs — that lie below the main sequence 
dwarfs on the HR diagram. 



The realm of the s ubdwarfs has been p reviou s ly exp lored bvlSandage fc EggenI (119591 ) . 
Hartwick et aP Jl984[ l. iMonet et all fll992h . loizisl Jl997l ). and I.Tao et al.l f l2n05h . to name a 
few. Subdwarfs' locations on the HR diagram are in part due to having metallicity abun- 
dances lower than most field stars, which causes their opacities to differ from those of reg- 
ular dwarfs. Subdwarfs are sometimes called low metallicity halo stars or Galactic thick 



disk stars based on their spectroscopic features, kinematics, and/or ages (IDigby et a. 



Lepine et al.ll2003l : ISurgasser et al.ll2003l : iReid and Gizid l2005l : iMonteiro et al.ll2006f ). Re- 



2003 



gardless of how they are described, subdwarfs are fundamentally different from their main 
sequence cousins. 

Because of their generally high intrinsic velocities in the Galaxy, many subdwarfs 
have been selected using high p roper rnotion efforts, such as the Lowell proper motion 
dGiclas et all Il97ll. Il978h. LH S Jhuvtenl ^9791 ). and Lepine-Shara Proper Motion- North 



(LSPM. iLepine fc Sharall2005al ) catalogs. Recently, subdwarfs have been selected by color s 
and spectroscopic observations via the Sloan Digital Sky Survey (SDSS, IWest et al.l 120041 ). 
After initial flagging as a potential subdwarf, spectroscopic and astrometric (i.e. trigonomet- 
ric parallax work) foUowup efforts are typically carried out to co nfirm or refute candidate s 
as true subdwarfs. Past subdwarf identi fi cation efforts i nclude iRvan and Norris 



Monet et all (119921). ICarnev et al.l (11994). iGizid (119971). iLepine et all (120031 ) . IWest et al. 



19911). 



kooi ). iReid and GizisI J2005h and iBurgasser fc KirkDatrickl tood). Most of these stud 



ies confirmed subdwarfs spectroscopically, but only lMonet et al.l (119921 ) provides the crucial 
trigonometric parallaxes that allow subdwarfs to be placed on the HR diagram. 



GizisI ( 119971 ) presented a pioneering effort to assign numeric subtypes for cool subdwarfs 



of spectral types K and M. First, he used the flux ratio of molecular band features, CaHn (n 
= 1-3; we will use "CaH" to indicate all three bands or index values throughout this paper, 
unless otherwise specified) and Ti05 with pseudo-continuum points to calculate spectro- 
scopic indices {hands/ f continuum)- On Ti05 vs CaHl or Ti05 vs CaH2-|-CaH3 plots, two 
high order polynomial lines can be used to separate a continuous distribution of stars into 
three categories: regular dwarfs, subdwarfs, and extreme subdwarfs. Numerical subclasses 
were then assigned using two independent linear fits for subdwarfs and extreme subdwarfs. 
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For a decade, this methodology has been used to assign spectral types for cool subdwarfs. 

However, the current method of assigning subdwarf spectral types is not directly linked 
to either their main sequence or giant counterparts, as is typically (but not always) the case 
with normal dwarfs and giants. In addition, we have found that many subdwarfs are assigned 
different subtypes even though their differences are limited to CaH. If these limitations can 
be overcome, a well-defined spectral sequence would benefit many research areas beyond 
classification efforts, including attempts to estimate effective temperatures and distances, as 
well as providing insight into understanding Galactic structure. 

In this work, we first provide clarification and recommendations for subdwarf termi- 
nology by addressing the usage of the confusing spectral prefix "sd". We then discuss our 
spectroscopic observations of 88 subdwarfs, generally targeting high proper motion stars in 
the southern sky. We next outline how synthetic spectra can assist us in understanding cool 
subdwarf spectral features. The bulk of this work describes a detailed effort to provide a 
subdwarf spectral sequence for stars having spectral types K3.0 through M6.0. We then 
apply this spectral typing method to those subdwarfs found in the SDSS database. Once 
we have a detailed understanding of what makes a star a subdwarf, we then discuss (1) 
why it is premature to assign precise parameters to subdwarf spectra, (2) misunderstandings 
related to the terms "extreme" and "ultra" subdwarfs, (3) why the old method works for 
dwarfs, but not sub dwarfs, and finally (4) the recent set of subdwarf spectral standards from 



Lepine et al.l (120071 ). 



'VI" Subdwarfs are Different from "sd" Subdwarfs 



The first subluminous objects fainter than main sequence stars were reporte d bvlAdams fc Jov 



19221) when they were trying to determine the luminosities of A- type stars. In lAdams et al. 



19351 ) they called these stars "intermediate white dwarfs" to separate them from typical 



white dwarfs, and in the same paper reported the first six "intermediate white dwarfs" — 
now known as LHS 405 (sdF3), LHS 540 (F8IV), LHS 1501 (A4p), LHS 2194 (sdF5), HD 
132475 (F5/F6V), and HIP 68321 (A4). All six have either A or F types (spectral types are 
from SIMBAD) in the modern MK spec t ral cl assification system. However, the term "sub- 
dwarfs" was not suggested until iKuipeii (119391 ). He expected three classes of objects to be 
found in his spectroscopic survey of high proper motion stars — white dwarfs, stars of large 
spectroscopic parallaxes, and a class that was 2 or 3 magnitudes less luminous than main 
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sequence stars of the same coloi0. He suggested the name "subdwarfs" be used to represent 
this final, independent, class of stars. This name paralleled the use of the term "subgiants" 
to describe stars that fall below the giants on the HR diagram. The name also eliminated 
the confusion with white dwarfs, which are much less luminous than main sequenc e stars 
(and we now know are a completely different type of object). A year later, iKuiped (119401) 
repor ted the first three M-type subdwarfs — Kapteyn's star, LHS 20, and LHS 64 (iGizis 
19971 confirmed that all are M subdwarfs). Although they were termed subdwarfs, he used 
the same spectral classification as dwarfs (MO and M2, see lKuiDei]ll940l Table 1). The "sd" 
spectral classification prefix for subdwarfs did not appear until I Joy Jl947h . when he used the 
strengthening of the Lindblad depression around 42 26 A. 

In the late 1940's through 60 's, the term "subdwarfs" also began being used for a 



class of underluminous blu e stars (IHumason &: ZwickyI 119471 : iFeigd Il958l : iGreensteinI Il966 



Greenstein fc Munchlll966l ). The terminology was based on the understanding that if these 
were high luminosity blue stars, their distances would be outside the Milky Way, so it 
was surmised that these stars should be underluminous, and therefore closer. Although 
their temperatures are similar to O and B dwarfs, their spectral features are, in fact, quite 
different. Generally, 0-type subdwarfs (sdO) and B-type subdwarfs (sdB) both have broad 
Balmer absorptio n. sdBs have weak o r no H e lines , while sdOs have strong He II (4686 A) or 
other He II lines (lHeberlll992l ). Since iFeigd (119581 ) . such blue objects have been called sdO, 
sdB, or sdOB-type stars. 

Thus, we are left with the unfortunate situation that there are two different classes 
of stars called "subdwarfs." One is located at the cool end of the HR diagram while the 
other is at the hot end. The two classes of stars are subluminous for completely different 
astrophysical reasons but sha re the same "sd" spectral classification prefix. Cool subdwarfs 
usually have low metallicity (jChamberlain fc AUeii Il95ll : iGreenstein fc Eggenl Il966l : iMould 
19761: lAUard fc Hauschildtl ll995l ). so their opacties are different from dwarfs. For example, 
AUard fc Hauschildtl Jl995h discussed the possible opacity sources for a solar type dwarf and 
a [m/H] =— 2.5 subdwarfs. TiO dominates the opacity sources in the optical band. However, 
because of a decreasing metallicity for subdwarfs, TiO opacity decreases dramatically. Hence, 
this less blanketing from TiO bands causes more continuum flux radiated from deeper and 
hotter layer of stellar atmosphere and their spectrum falls closer overall to that of a black- 
body, so these subdwarfs appear bluer than dwarfs, as shown in Figure [^. To the contrary. 



^Kuiper assumed all of these high proper motion stars had Vtan less than 474 km/sec. Under this 
assumption, he could then make crude estimates of their absolute magnitudes without having trigonometric 
parallaxes. The results placed some stars three magnitudes below the main sequence. 

^The details of these synthetic spectra, the "GAIA model grids", are discussed in section 4. 
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hot blue subdwarfs of types O and B are progenitors of white dwarfs, and their subluminous 
nature is not caused by metaUicity at all. Instead, hot subdwarfs represent a stage in the 
stellar evolution cycle of an evolved star, and they happen to be crossing the main sequence 
at the moment of observation. 

Because they are different kinds of stellar objects, we su ggest that the two classes 
should not share the same spectral classification notation, "sd." iRomanI (119551 ) argued that 
for types later than GO, the spectral notation "VI" should be u sed for stars that a re 



2 magnitude less luminous than main sequence stars. Although IJaschek fc Jaschekl (119871 ) 
stated that "this designation (VI) should definitely be abandoned", no specific reasons were 
actually given. Here we propose that the luminosity class "VI" should be adopted for cool 
subdwarfs, especially for K and M types. The three primary reasons are: 



Cool subdwarfs of types K and M lie clearly below main sequence stars on the HR 
diagram, forming an additional class of objects. Assigning their types as VI continues 
the progression outlined by other classes. Like the family of giants, which includes 
luminous supergiants (I), bright giants (II), normal giants (HI), and subgiants (IV), 
the family of dwarfs includes main sequence dwarfs (V) and their subdwarf (VI) coun- 
terparts. 

Cool subdwarfs and OB subdwarfs are completely differently types of objects with 
different origins. The notation "sd" is suitable for OB subdwarfs, which have not yet 
as a group been established to have any sort of sequence on the HR diagram, at least 
when using using parallaxes from the Yale Parallax Catalog, Hipparcos, and more 
recent efforts. On the other hand, there are dozens of cool subdwarfs with parallax 
measurements, and they do form a coherent group below the main sequence on the HR 
diagram, making a Roman numeral designation reasonable. 

Historically, the Roman numerals used for luminosity categorization track with dif- 
ferent gravities for the stars classified, with higher gravities being assigned higher 
Roman numerals. We find that assigning "VI" for subdwarfs appropriately contin- 
ues this trend. Figur e [2] shows the mass-gravity relation using data from Table 1 of 



Lopez-MoralesI (120071 ). In the best represented mass regime from O.35M0 to O.7OM0, 
a crude trend indicates that lower metaUicity stars do indeed have higher gravities. 
Although much more data are needed to understand clearly how metaUicity affects the 
mass-gravity relation, current evidence supports using the "VI" designation for low 
metaUicity subdwarfs, which tend to exhibit higher gravities. 



Thus, in order to separate the OB subdwarfs from the cool subdwarfs, we suggest the 'sd" 
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prefix should not be used for low metallicity (and/or high gravity, as will be shown below) 
sub dwarfs. 



3. Observations and Reductions 



3.1. Observations 



Our subdwarf targets were sel ected from several different efforts , including lists of spec- 
troscopically identi fied subdwarfs ( Gizis Il997 : Reid and Gizis 2005 ). subdwarfs with paral- 



lax measurements (Jao et al 



catalog s (IDeacon et al. 



20051: 



2005 



Costa et a 



Subasavage et al 



12005 



2005a 



20061 ). candidates from proper motion 



b|), and stars with metallic ity measure- 



ments ( ICarney et al.lll994l : ICayrel de Strobel et al.ll200ll : iNordstrom et al.ll2004l ). We define 
our subdwarfs of interest to be those with V — ^ 2.0, [m/H] ^ —0.5, or having absolute 
magnitudes at least one magnitude less luminous in Mks than a fit to main sequence stars 
of comparable color with trigono metric parallaxes from the RECONS (Research Consortium 
on Nearby Stars) 10 pc sample (IHenry et al.l 120061 ). 



Spectroscopic observations were made with the 1.5-m and 4.0-m telescopes at CTIO. 
For the observations on the 1.5-m from 2002 to 2006, the R-C spectrograph with a Loral 
1200x800 CCD camera was used with the #32 grating (in first order) at tilt 15.1°. The 
order-blocking filter OG570 was utilized to provide spectra covering the range of 6000A to 
9500A with a resolution of 8.6A. The only variation in observing parameters was that during 
the May and December 2006 observing runs a larger slit width of 6" and 4", respectively, 
was used instead of the 2" slit used in previous runs, in order to minimize the differential 
color refraction (because the slit orientation was not changed during observations). For 
observations on the 4.0-m in 2002, the R-C spectrograph with a Loral 3KxlK CCD was 
used with the #181 grating (in first order) at tilt 58.8°. The order-blocking filter OG515 was 
utilized to provide spectra covering the range from 5500A to lOOOOA with a resolution of 6A. 
Fringing at wavelengths longer than ~7000A in the 4.0-m data was removed by customized 
IDL routines. Bias frames and dome fiats (and sky flats at the 1.5-m) were taken at the 
beginning of each night for calibration. At least two exposures were taken for each object 
to permit cosmic ray rejection. If stars were faint, additional observations were sometimes 
made. A 10 second Ne+He+Ar or Ne only arc lamp spectrum was recorded after each 
target to permit wavelength calibration. Several spectroscopic flux standard stars found in 
the I RAF spectroscopy reduction packages were observed during each observing run, usually 
nightly. Reductions were carried out in the standard way using IRAF reduction packages. 
Wavelength and flux calibrations were done using onedspec.dispcor and onedspec. calibrate 
within IRAF, respectively. 
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Many of the subdwarfs discussed in this paper have new trigonometric parallaxes and 
VRI photometry acquired during our southern nearby star program, C TIOPI (Cerro Tololo 
Inter-american Observatory Parallax Investigation, see lJao et al.ll2005l ). In a future paper in 
this subdwarf series, we will present the astrometric (particularly trigonometric parallaxes) 
and photometric results. 



3.2. Identifying Subdwarfs 



During our five year spectroscopic campaign, we have acquired spectra for more than 
900 objects. To glean subdwarfs from our spectroscopic database, we calculated their spec- 
troscopic in dices l i sted i n Tabled] (targets listed alphabetically), mimicking the methodology 
outlined bv lcizisl Jl997h . 



Figure [3] shows Ti05 plotted against CaHl and CaH2+CaH3 for various samples of 
stars. Our subdwarfs are shown with solid circles. For comparison, small dots indicate 



main sequence stars from iHawley et al.l (I19961). while o pen t riangles and squares r epresent 
subdwarfs and "extreme" subdwarfs from iGizid (119971 ) and iReid and Gizij ( 120051 ). Some 
stars (solid circles) having indices located near or in the main sequence regions in these plots 
have been manually checked to confirm that they are subdwarfs. Using the HR diagram 
in Figure HI we confirm the low luminosities of our spectroscopically selected subdwarfs 
{V — Ks > 2.7) that have accurate trigonometric parallaxes. 

Although we focus primarily on the K and M subdwarfs for this study, we include a 
few G-type subdwarfs among those selected from metallicity measurements in the literature. 
For reference, we consider early K-type stars to have types KO to K2, mid K-types to be 
K3 to K5, and late K-types to be K6 and later. We find that it is difficult, but possible, to 
separate late G from early K-type stars using our spectral coverage and resolution. Their 
continuum slopes have only slight differences across our wavelength window coverage, and 
there are no notice a ble ab sorption differences bey ond 7500A. Spectra for GIV to K5V types 
from Ijacoby et al.l ( 19841 ) (resolution ~4A) and ISilva fc Cornelll ( 19921 ) (resolution ~llA_) 
are plotted in Figure [51 The only strong features are the absorption lines of Ba I (6497A)o 
and Ha (6563A), with gradually increasing Ba I absorption and decreasing Ha absorption 
as the effective temperature drops. These effects can be seen in both sequences, regardless 
of the spectral resolution. Our spectral resolution of 6-9A falls between the resolution of 



Turnshek etHI |l985h noted that this feature at 6497A is a blend of different atomic lines, including Fe 



I, Ba I, Ca I, Mn I, Co I, Ti I and II, and Ni I. Ba I is likely the dominant absorber because it has the largest 
Einstein coefficient in the NIST atomic spectra database. 
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the two sequences shown, so we can use the relative absorption strengths of Ba I and Ha to 
separate G and K-type stars. In total, we have identified 88 K and M subdwarfs and five G 
type subdwarfs, using spectra with coverage from 6OOOA to 9OOOA. 



3.3. Sorting Spectra 

After reduction, the spectra were sorted into different bins based upon similarity in 
overall slope and features. This assured that stars in each bin had approximately the same 
temperature. However, several impediments to clean sorting were encountered: 



• All mid K-type subdwarfs had spectra virtually indistinguishable from dwarf standard 
stars, yet they had low metallicity measurements and/or were found below the main 
sequence on the HR diagram (see top figure in Figure [6]). 

• Many subdwarf spectra placed into the same bin showed differences only in CaH (see 
middle figure in Figure [6]) . 

• Many subdwarf spectra matched different dwarf standards at the blue (A < 7570A) 
and red (A > 8200A) ends (see bottom figure in Figure [6]). 



In order to understand what factors caused these anomalies, we next examine theoretical 
studies that provide synthetic spectra that can be compared directly to the observed spectra. 



4. Grids of Synthetic Spectra 



We use grids of synthetic spectra computed with PHOENIX codes (hereafter, GAIA 

model grids) to understand how subdwarfs' physical parameters (temperatu re, metallicity 

and g ravity) affect their spectra. The most recentl y relea s ed GAIA model grids (IBrott fc Hauschildt 
2OO5I) are available at an FTP site i n HamburJ^- iGizisI fll997h . IWoolf fc WallersteinI j2o"o5f L 



and iBurgasser fc Kirkpatrickl (120061 ) have all used these synthetic model grids to character- 
ize subdwarfs, but an older version of the grids was used in all three cases. The version we 
employ here, 2.6.1, was released in late 2004. A comparison of one pre-2004 spectral model 
(provided by V. Woolf, private communication) and a new spectral model (from the GAIA 
model grids) for a cool subdwarf is shown in Figure [71 Improvements to the new models 



''ftp://ftp.hs.uni-hamburg.dc/pub/outgoing/phoenix/GAIA/v2.6.l/ 
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include (1) an enlarged and enhanced version of the equation of state, (2) more atomic, ionic, 
and molecular line opacities, (3) inclusion of the formation of dust particles for cool stars, 
and (4) microturbulence calculations. Additional water and TiO opacities and the inclusion 
of dust are enhancements particularly applicable to the low mass stars discussed here. 

The two spectra shown in Figure [7] are virtually identical redward of 7000A, but there 
are significant differences between 6500A and 7000A, where the CaH2 and CaH3 bands are 
found, and these differences will certainly affect evaluations done with the older models. The 
new model in this region has much shallower absorptions than the old one, which will affect 
metallicity and gravity estimates. A few narrow absorption features (Li I at 6103A, Ca I 
at 6122A and 6162A) are also changed. As outlined in the discussion section ^ even the 
latest version of the model grids does not provide ideal matches to real spectra, so further 
progress can still be made. 



4.1. Identifying Mid K-type Subdwarfs 

We use the latest GAIA model grid^ to calculate predicted CaH band strengths and 
plot the derived indices against effective temperatures in Figure [HI We evaluate stars with 
Teff between 2700K and 4500K and [m/H] between 0.0 and —3.0. For the moment, we 
adopt log g = 5.0 generically for subdwarfs. This gravity value does not apply to all types 
of subdwarfs, but we are presently interested in outlining the behavior of the CaH features 
with metallicity alone. 

When Teff is less than about 3300K, Figure M shows that the CaHl index decreases 
(stronger absorption) when metallicity decreases from 0.0 to —2.0 at fixed Teff. However, 
the trend reverses for [m/H] = —2.5 and —3.0. When Teff is between 3300K and 3500K, 
there is a very weak relation between metallicity and the CaHl index, but from 0.0 to —1.0, 
the relation (CaHl index decreases when metallicity decreases) still holds. This relationship 
is degenerate for lower metallicities. For temperatures hotter than 3500K, the CaHl index 
increases (weaker absorption) as metallicity decreases, in contrast to the low temperature 
region. For the CaH2+CaH3 index, the trends are generally the same, except that (1) 
at temperatures less than 3200K the index decreases (stronger absorption) only for 0.0 to 
— 1.0, with a reversal for lower metallicities, and (2) the trend for higher temperature stars 
(increasing index, weaker absorption with lower metallicity) is the same, but the turnover is 



^GAIA model grids also provide various valu es of g-elements (O , Ne, Mg, Si, S, Ar, Ca, and Ti), which 
yield abundance ratios such as 0/Fe, Ne/Fe, etc. iNissen fc Schusteii (jl997T ) showed that most field halo stars 
have [a/Fe]*«[Q;/Fe]0, so we select [Q;/a0]=O.O. 
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near 3200K rather than 3500K. 

Even more important than these subtleties is that overall, the hotter the subdwarf, the 
weaker its CaH bands. Note the collapse of any differences between indices for hotter stars 
of various metallicities in Figure [H This collapse makes separating mid K-type subdwarfs 
from dwarfs based only on spectroscopic indices difficult using our spectral coverage (6000A- 
9000A) and resolution (6A-8.6A). An alternative method, such as the HR diagram shown 
in Figure HI sufficiently solves the problem for mid K-type subdwarfs and will be discussed 



Bonfils et al. 


2005; 


Bean et al. 


2006) 



not of sufficiently high resolution to measure metallicities, but all selected K-type subdwarfs 
have measured [m/H] ^ —0.5 from other publications. 



4.2. Mid K-type Subdwarf Sample from Our Observations 

There are 31 subdwarfs having CaHl indices larger than 0.9 enclosed by the dashed 
box in Figure O In this region, there are no subdwarfs with previously measured CaH/Ti05 
indices to compare to our new sample of mid K-type stars, making it difficult to separate the 
dwarfs and subdwarfs based on the CaHl index alone. Fortunately, 21 of these 31 stars have 
trigonometric parallaxes, and are plotted with concentric circles in Figure HI At least 15 of 
these stars are subdwarfs based on their locations one or more magnitudes below the main 
sequence on the HR diagram. The star with the largest offset is DEN0515-72lfl located at 
{V — K) = 3.3, Mks = 9.7. This star is a full 4.5 magnitudes less luminous than the main 
sequence, but has no CaH or Ti05 features. Two additional stars (G016-009AB and G026- 
009ACD) above the main sequence are known to be double-line spectroscopic binaries with 
[m/H] < —0.5. These individual targets are discussed in section 6.2. The dwarf/ sub dwarf 
status of only four stars of the 21 remain ambiguous — we suspect that most of them are also 
subdwarfs, perhaps with as yet undetected companions elevating them into main sequence 
territory. 

Generally, the spectroscopic index method fails to distinguish subdwarfs from dwarfs if 
the derived CaHl index is greater than 0.9. With the benefit of additional trigonometric 
parallax information and/or metallicity measurements, however, we can conclude that nearly 
all of these mid K-type stars are indeed subdwarfs. 



^This star was first reported in ICosta et al.l (2006!) reference star #4 in the LHS 1749 parallax field. 
We identify it henceforth as DEN0515-7211 (DENIS-P J051545. 1-721122). 
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5. Late K-type and M-type Subdwarfs from GAIA Model Grids 

For late K-type (redder than K5.0) and M-type stars, we use the GAIA model grids 
to understand how the effective temperatures, metallicities, and gravities affect the shapes 
and features of subdwarfs' spectra. This analysis allows us to develop a spectral sequence 
for cool subdwarfs of types K6.0 to M6.0, which is tentatively extended blueward to K3.0 
when additional information is incorporated. We plot GAIA synthetic noiseless spectra in 
Figure [H at increments of 200K (cooler than 4000K) and 400K (hotter than 4000K) for stars 
with metallicities of [m/H]= 0.0, —1.0, and —2.0. Because gravities have very limited impact 
on the overall shapes of the spectra (shown in the top panel of Figure [TUi) . we do not show 
gravity plots with fixed temperatures and metallicities (of course, some features do change 
markedly with gravity, but not the overall slopes of the spectra). Based on the synthetic 
spectra from GAIA model grids, cool subdwarf spectra between 6OOOA and 9000A exhibit 
the following trends: 



The effects of metallicity are minimal in low resolution subdwarf spectra for stars with 
temperatures of 4400K and hotter. This makes it difficult to separate dwarfs and 
subdwarfs using low resolution spectra (as discussed in §4.11) . However, from these 
noiseless spectra, we can still identify a few metallic lines showing metallicity trends. 
The most prominent feature is marked #1 at 6256 A in Figure [H Unfortunately, we 
do not see this featu re in any of our K-type spectra (nor is it listed in Table 1 of 



Turnshek et al.lll985l ). so we consider its validity questionable. The next prominent 
feature, marked #2 in Figure [HI is Ca I (6162A). This line can possibly be used to 
distinguish subdwarfs from dwarfs (see examples discussed in §7.2p . but in practice it 
is somewhat difficult to evaluate in real spectra (with noise) at our resolution. 

For stars with temperatures of 2800-4000K, metallicity strongly affects the spectra 
between 6OOOA and 8200A. This is the region that has been historically used to assign 
spectral types. In effect, subdwarfs with decreased metallicities have spectra that are 
"brightened" or "less blanketed" at the blue end, relative to solar metallicity stars. 
However, the continuum at wavelengths longer than 8200A for temperatures 3400- 
4800K is nearly free of metallicity effects. We can therefore use the 8200-9000A region 
to establish subtypes in the subdwarf spectral sequence because the slope is a function 
of temperature. This also allows us to mirror the spectral sequence for dwarfs, providing 
a useful link between the dwarf and subdwarf sequences. 

The Ti05 band strength at 7050-7150A is very sensitive to metallicity for temperatures 
cooler than 4000K. As shown in the top panel of Figure [TOl the Ti05 band strength is 
effectively independent of gravity. We can therefore use the Ti05 feature to separate 
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subdwarfs with different metallicities if their continua (8200-9000A) are the same, 
regardless of their gravities. 

• For a star of given temperature, stronger CaH bands could be caused by lower metal- 
licity, as shown in Figure [9l or higher gravity, as shown in the top panel of Figure [TOl 
If two subdwarfs have the same continua from 8200-9000A and the same Ti05 band 
strength, but their CaH bands are different, we can rank them by their relative gravi- 
ties. 

Consequently, the impediments to sorting cool subdwarf spectra at our resolution dis- 
cussed in §3.31 can be overcome by understanding the trends revealed in GAIA synthetic 
spectra: (1) the mid K-type dwarfs have the same spectra as dwarfs for our spectral cover- 
age and resolution, (2) CaH features are affected by both metallicity and gravity, while the 
Ti05 band is affected by metallicity but not gravity, and (3) the continuum from 8200A to 
9000A is not strongly affected by either metallicity or gravity, so can therefore be used for 
spectral sequencing. 



6. Additional Evidence Supporting the Metallicity and Gravity Trends 

Indicated by the GAIA Models 

To investigate the metallicity trends seen in GAIA models, we compare our available 
spectra to metallicity measurements provided independently by others. Measuring M dwarf 



metallicities is difficult, bu t several recent attempts have made pr ogress (IValenti et al.lll998 



Woolf fc Wallerstein 2005; Bonfils et al. 2005 



Bean et a. 



have observed are included in the study by iBonfils et al 



20061) . Six red dwarfs that we 
( 20051 ). comprising three pairs of 



M dwarfs of types Ml.OV, M2.5V and M3.0V shown in Figure [TTl Each pair includes a 
relatively low metallicity dwarf (gray), and a relatively high metallicity dwarf (black). The 
red ends (from 8200A to 9000A) of each pair match one another, but the blue ends of the 
lower metallicity members is brighter in each predicted by the GAIA models. The 

effect is rather subtle, but the derived metallicities for each pair are not wildly different (none 
of the six stars are subdwarfs), and yet the trend is confirmed in all three cases. One caveat 
is that the metallicitie s from these six stars were determined from the polynomial relation 



m 



Bonfils et al.l (120051 ) . not measured directly from spectra, but this appears to be the best 



that can be done given the available data. We conclude that the metallicity trend revealed 
in the GAIA models (see Figure ^ is sound because it appears to be confirmed in real M 
dwarf spectra. 



Contrary to the evidence for the metallicity trend, we have found no direct spectroscopic 
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results to support the gravity trend in dwarfs and subdwarfs. Direct gravity measurements 
are difficult because (1) stars must have both mass and radius measurements and (2) clean 



spectra without contamination from companions must be obtained. This 



imits the available 



target hsts to eclipsing binaries such as those discussed in iLopez- Morales! (120071 ) with cleanly 
deconvolved spectra or visual binaries in which individual radii can be measured via long 
baseline interferometry. Neither class of objects yet provides a rich dataset for cool dwarfs 
or subdwarfs. 

Hence, we rely on other observational or theoretical efforts to investigate the gravity 
trend. The top panel of Figure [10] indicates that other than some sharp metallic lines, CaH 
bands sh ow the most pr ominent changes when gravity varies (TiO is unaffected by changing 



gravity). lOhmaru (119341 ) demonstrated that the CaH2 band is found in the spectra of M-type 



dwarfs, but is not observed in the spectra of M-type giants, thus identifying the CaH 2 band 



as a g ravity indicator to separate dwarfs and giants. In addition to the GAIA models, iMould 



(119761 ) also showed that for stars with effective temperatures of 3250K, a spectrum from his 
atmospheric model with log g = 5.75 has stronger CaH2 than a spectrum with log g = 4.75. 
If the CaH2 band is a gravity indicator, we may presume that the same gravity effects for 
CaHl and CaII3 bands will be seen. 

As shown in the two panels of Figure [10], the GAIA models imply that the CaH band 
strengths are indicators of both gravity and metallicity differences. In reality, if two red 
subdwarfs have spectra with the same overall continua and slopes, as well as matching TiO 
band strengths (e.g. middle panel in Figure [H]), the only remaining discrepancies will be at 
the CaH bands. We believe that such differences are caused by different gravities. 

Obviously, there is not yet a wealth of accurate direct observational results of red dwarfs 
and subdwarfs that can be used to stress test the metallicity and gravity trends seen in GAIA 
models. Nonetheless, what little we do have supports the trends, so we use these trends to 
assist us in establishing the subdwarf spectral sequence discussed below. 



7. Subdwarf Spectral Sequence 
7.1. Procedures 

In addition to presenting 88 subdwarfs, a goal of this project is to establish a subdwarf 
spectral sequence that mirrors the sequence for dwarfs and considers the latest available 
synthetic models. Although synthetic models are not yet capable of fully representing the 
complicated spectra of these cool stars, the models can be used to investigate the primary 
factors that affect subdwarfs' spectra. Specifically, we examine the spectra in the frame- 



-14- 



work of what appear to be the three main drivers of the trends observed — temperature, 
metalhcity, and gravity. 

We first separate our available sub dwarf spectra into several groups that have similar 
overall slopes. Known subdwarfs such as LHS 12, GJ 161, and LHS 2734A are used as anchor 
points. Within each group, several stars that have high surface gravities are obvious because 
their spectra match except in CaH, implying nearly identical temperatures and metallicities. 

We then compare each subdwarfs spectrum with our sequence of dwarf spectral stan- 
dard^ that span types from KO.O to M9.0 using an IDL program to find the closest match 
to the continuum slope in the region 8200-9000A. Visual checks of the matches between all 
subdwarf and standard dwarf spectra are also made to ensure match quality. 



7.2. Results and Notes on Objects 

Based on the 88 confirmed cool subdwarf spectra we have, we present a sequence of 
subdwarfs with spectral types spanning K3.0 to M6.0, listed in Table [21 We also identify five 
additional G type subdwarfs. After the five G-type subdwarfs, we sort the cool subdwarfs 
from K3.0[VI] to M6.0VI, using double lines in Table [2] to separate each type. 

As discussed in §4.11 because mid K subdwarfs are virtually indistinguishable from K 
dwarfs at our spectral resolution, we use [VI] to indicate their questionable luminosity classes, 
which are currently based on their metallicities, kinematics, or locations on the HR diagram. 
We anticipate that higher resolution spectra will reveal these stars to be subdwarfs. A colon 
after the type indicates that we have had difficulty in assigning a subtype, metallicity, or 
gravity. 



^ Our dwar f stand ard sequence is a hybr i d of c ool dwarf standards from Gray et al. ( 2006h for K stars, 



and Boeshaar ( 1976h and Kirkpatric 



a future publication (|Beaulieu et al.l 
sequence from SDSS 



c et al.l (|l99lh for M stars. The dwarf standard sequence is the topic of 
200 9|). Our spectral sequencing efforts began before the recent spectral 



Bochanski et al was released. We continue to use our standard sequence because 



(1) stars in our sequence have trigonometric parallaxes so we can understand how metallicities and gravities 
affect stars' positions on the HR diagram, (2) the SDSS sequence do es no t include K - type stars, (3) we 
have acquired spectra for many K-type stars from Gray et al. ( 20031 ) and Gray et al. ( 2006 ) so we have 
benchmark K dwarf spectral standards to link to the M dwarf sequence, (4) SDSS standard stars have telluric 
lines removed (It is difficult for us to make comparisons to SDSS spectra in the 6800-7100A region, which 
includes O2 and H2O absorption. This region overlaps with the CaH3 and Ti05 bands. Thus, all subdwarfs 
would appear to have high gravity when compared to the SDSS standards.), and (5) our dwarf spectra have 
been acquired using the same telescope/instrument/observing protocols as used for the subdwarfs, and the 
data have been reduced identically, so they are systematically consistent. 
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Within each type for which sufficient spectra are available, we sort targets by their 
metallicities (lowest metallicity ffist). A letter "m" indicates a star having the same metal- 
licity as a main sequence dwarf, while more negative signs indicate lower metallicities, e.g. 

m is more metal poor than m— . We use as many as six negative signs, because in the 

case of the Ml.OVI type, we have 23 stars that fall into seven different metallicity categories. 
Gravities are indicated by "g" with additional plus signs for higher gravities, e.g. g++ is 
higher gravity than g+. A baseline subdwarf of low metallicity is assigned m— and g. If a 
subdwarf has a similar metallicity but higher gravity, the designations are m— and g+. A 
few stars appear to have solar metallicity and are subluminous only because of high gravities; 
these stars are assigned m and g+. We have removed any stars from this study that might 
have gravities lower than main sequence stars, i.e. slightly evolved stars such as subgiants, 
that would have gravity g— . 

Note that values for metallicity and gravity are not comparable across all spectral types, 
i.e. Ml.OVI with metallicity "m— " is not equivalent to M2.0VI with metallicity "m— ," nor 
is "g+" for Ml.OVI the same as "g+" for M2.0VI. We can hope to formalize a definitive 
subdwarf spectral sequence that includes temperature, metallicity, and gravity trends when 
hundreds of systematically consistent subdwarf spectra and improved model grids are avail- 
able. Objects with virtually identical spectra are listed in alphabetical and/or numerical 
order. 

We discuss each spectral subtype from K3.0[VI] to M6.0VI in the following sections and 
highlight noteworthy subdwarfs of various types. The order of the highlighted targets is 
based on their metallicities or gravities. Colored spectra in the Figures in these sections are 
illustrative only — they do not represent any numeric metallicities or gravities, so a red-lined 
spectrum in one type does not necessarily have similar attributes to a red-lined spectrum in 
another type. 



7.2.1. K3. 0[ VI] to K5. 0[ VI] types 



We begin with K3.0[VI], for which it is still difficult to separate subdwarfs from dwarfs 
at our spectral resolution. Figure [9] shows how similar mid K-type subdwarfs and dwarfs 
are (temperatures 4400K and 4800K). For these stars, we use other independent published 
measurements, e.g. [m/H] values, to con&m their subdwarf natures. We currently call stars 
with [m/H] < —0.5 subdwarfs. The K3.0[VI] to K5.0[VI] spectra are shown in Figure] 



G 016-009 AB (K3.0[VI]) IColdberg et all (120021 ) reported this star to be a double- 
lined spectroscopic binary with P=9.9 days, T^ff = 5903K, and having a mass ratio (M1/M2) 
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of ~1.22. Given that the components are presumably coeval, the combined spectrum should 
represent the metallicity for each component. Although their estimated temperature i s that 
of a G type star, our spectr um has the same slope as a K3.0V. ICayrel de Strobel et al.l (I2OOII ) 



and iGoldberg et a 



(I20021) reported measurements of [Fe/H]= —0.7 and —1.0 respectively, 

while Ihaird et al.l JlQSsh found [m/H]= —1.11. Apparently, it is a low metallicity star and 



its location on the HR diagram is elevated becaus e of multiplicity. This system is also 
reported to be a photometrically variable system in iKazarovets et al.l (120061 ). This system 
demonstrates that low metallicity subdwarfs do not have different spectra from their main 
sequence counterparts at our wavelength coverage and resolution, so we assign its luminosity 
class as [VI] as a result its low metallicity measurements. 

LHS 2467 (K 4.0[VI1) Our spe ctrum shows that it is a K4.0[VI], rather than a G7V 
type, as reported in iBidelmanI (119851 ). Its continuum is not as blue as KO.OV and its Ba 
I/Ha lines are not a s weak /stro ng as shown in Figure [5] f or a G type star. Its weighted mean 
parallax from lESAl (119971) and Ivan Altena et al.l (119951 ) is 10.63± 1.88 mas and its proper 



motion is 0'.'98 yr~^ ( lLuytenlll979l ). indicating Vtan = 437 km sec~^. iRyan and Norris I (Il99ll ) 
reported the star to have Vrad = 4:4: km sec~^. After removing the solar motion, LHS 2467 
has {U,V,W) = (281.4,-304.8,-41.9) km sec~^. The extremely high tangential and space 
velocities are indicative of a subdwarf. 

G 026-009ACD (K5.0[VI]) The wide (132", corresponding to 6520 AU at the sys- 
tem's distance of 49.4 pc) common proper motion companion, B, in this quadruple system is 
a white dwarf known as G026— 010 or G026— 009B, that during its planetar y nebula phase 



should have had limited impact on the metallicity of the close ACD triple. iPeterson et al. 



( 1l980l ) reported G026-009AC t o be double- l ined s pectroscopic binary with P=3.75 days, with 
a mass ratio (M^ / My)^1.25. Allen et al. ( 2000) reported a third component, D, 0'.'7 away 



from G026-009AC. iMorrison et all (120031 ) and lAllen et al.l (120001 1 both report the system 
(ACD) to have low metallicity, measuring [Fe/H]= —0.91 and —1.19, respectively. Both the 
V and Kg magnitudes include three stars, causing the point to be improperly placed on the 
HR diagram. 



G 022-015 (K5.0[VI];) IWoolf fc WallersteinI (120051 ) reported this system to have [Fe/H] = 
—0.61. G022-015's spectrum is virtually identical to G026-009ACD, as shown by plotting 
both together in Figure | 



GJ 223.1 (K5.0[VI]) This star has the same slope as the K5. 0V standard and its 
location on the HR diagram is on the main sequence line. Nonetheless, IWoolf fc Wallerstein 
( I2OO5I ) reported it to have [m/H]= —0.62, so we tentatively consider it to be a subdwarf 
until further information indicates otherwise. 
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Although not definitive, the Ca I line at 6162A in these spectra seems to show the 
metalhcity trend seen in the noiseless model spectra plotted in Figure [9] and discussed in 
section O LHS 2467 has a much weaker Ca I line, and G016-009AB and G026-009AB have 
slightly weaker Ca I lines than seen in the dwarfs' spectra. However, G022— 015 and GJ223.1 
both have approximately the same Ca I line strengths as dwarfs. Contrary to the Ca I line, 
the Ca II lines at (8542A and 8662A) do not show any metallicity trend from our spectra. 
Because these prominent lines fail to separate subdwarfs from dwarfs at our resolution and 
S/N, other evidence, i.e. kinematics, parallaxes or independent metallicity measurements, 
are required to identify the subdwarfs. An alternative method to separate the subdwarfs 
from the d warfs is to ob tain spectra that include the MgH bands at 4845A , 521 lA, and 



562lA that iBesselll ( 19821 ) pointed are sensitive to metal abundance for temperatures hotter 
than 4500K. 

Several possibilities may explain the locations of the three subdwarfs, LHS 2467, G 
022-015, and GJ 223.1 on the HR diagram just below or on the main sequence line: (1) they 
are subdwarfs with unseen companions that brighten their Mks magnitudes, (2) they are 
slightly evolved and have moved significantly from the subdwarf zero age main sequence line, 
or (3) their metallicity measurements are not accurate and they are, in fact, main sequence 
dwarfs. 

LHS 125, LHS 232 and LHS 327 The spectra for these three stars, shown in Fig- 
ure [121 match neither the LHS 2467 K4.0[VI] spectrum nor a K4.0V spectrum, although the 
spectra have the overall slopes and the stars have V — Kg colors indicative of K4.0 stars. 
They are tentatively assigned types of K4.0[VI] based on their locations on the HR diagram. 



7.2.2. K6.0VI 

We identify four stars to have spectral type K6.0VI, shown in Figure [HI Their spectral 
slopes fall between K5.0V and K7.0V. 

LHS 193 A This is a binary with a separation of 12'.'59 co mprised of a cool subd warf 



(A) and a DC- type white dwarf (B) with a featureless spectrum (IMonteiro et al.ll2006l ). At 
the system's distance of 31.2 pc, the large separation of the pair corresponds to 393 AU, 
indicating that significant pollution of the subdwarf by the white dwarf during its planetary 
nebula phase seems unlikely. From Figure [TU it is clear that LHS 193A's spectral slope 
is between K5.0V and K7.0V. Although the spectrum is a near match to K5.0V between 
8200A and 9000A, the blue end (6000A-7500A) is too low to be a K5.0V (as a subdwarf, 
the blue end of the spectrum would have to be above the dwarf standard). As shown in 
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Figure m for a star with 4400K (a K5.0V from lAUen's Astrophysical Quantitiesll2000l ). there 
is effectively no difference in the continua of spectra for dwarfs and subdwarfs. Hence, we 
assign LHS 193A a type of K6.0VI because its Ti05 band is "brightened" (low metallicity), 
as shown in Figure M Given its low metallici ty, location on t he HR diagram (shown in the 
inset of Figure [H]), its Vtan =147.3 km sec~^ Jjao et al. 2005), and its age of 6-9 Gyr based 
on its white dwarf companion's cooling age (IMonteiro et al.ll2006l ). it is most likely a thick 
disk sub dwarf. 

LHS 73 and LHS 227 These two objects have virtually identical spectra, so we only 
show LHS 227 in Figure [H] for clarity. The only differences between these spectra and that 
of LHS 193A are in CaH. Using results from the top panel of Figure [TOl we find that LHS 
73 and LHS 227 have higher gravities than LHS 193A. 

We have acquired spectra for both LHS 72 and LHS 73, a wide (~96" separation) 
common proper motion pair. Unfortunately, the spectrum of LHS 72 {V — Ks = 3.29, 
Mks = 6.69) is poor and requires re-observation, although it would fall in the subdwarf 
region of the HR diagram. LHS 73 has V — K^, = 3.42 and Mj^,^ = 7 .27 and is plotted in 



Figure UM also clearly in the subdwarf region. iRodgers fc Eggenl (Il974r) previously identified 
both objects as subdwarfs, but no spectral types were given. iBidelnianI (|l985l) identified 
LHS 72 and LHS 73 as K4 and K5 dwarfs, respectively. iReyle et al.l (120061 ) reported LHS 
72/73 (they identified LHS 73 as sdK7) to be the nearest subdwarf binary based on their 
spectroscopic parallax (18.7 pc), although the trigonometric parallax from YPC is 37.6±8.9 
mas (21.5 pc < d < 34.8 pc). The nearest known s ubdwarf bin ary system is actually /x Gas 
AB, at a distance of 7.5 pc {ntrig = 132.4±0.6 mas. IeSA|E997[ ). 



LHS 161 This star was previously reported to be an esdM2.0 extreme subdwarf (IGizis 



19971 ). However, its slope and spectral features, like LHS 227, are too hot for type M2.0VI. 
As shown in Figure[T3l LHS 161 has stronger GaH than LHS 193A and LHS 227, but the rest 
of the spectrum is the same. We conclude that LHS 161 has higher gravity than the other 
two subdwarfs. LHS 161 's very strong GaH lines and consequent GaH and Ti05 indices 
place it in the previously called "extreme" subdwarf region, which typically implies that the 
star has very low metallicity. It appears, from our current understanding, that the strong 
GaH in LHS 161 may not be linked to metallicity, because the Ti05 band is not significantly 
different from LHS 193A, LHS 73 or LHS 227. 

In the inset of Figure [H we connect LHS 193A, LHS 73, LHS 227, and LHS 161 on the 
HR diagram with thick arrows to outline a sequence of increasing gravity effects that shift 
their locations to the lower right (less luminous and redder). We will soon see more such 
gravity effects, which are always towards lower luminosities. 
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1.2.3. K7.0VI 

We identify four stars to have spectral type K7.0VI, shown in Figure [151 This is the 
earhest spectral type for which we see clear metallicity effects in the blue regions of our 
spectra. 

LHS 2734A This is the primary in a new common proper motion system discovered 
during CTIOPI, for which we measure the components to have /i = 0'.'59 and 0'/60 yr~^ 
at position angles 228.2° and 228.4° for the A and B components, respectively. The B 
component is located 68'.'8 away at position angle 162.4°. Unfortunately, both components 
are beyond our current trigonometric parallax limit (100 pc), but the A component provides 
a useful subdwarf anchor point even though its location on the HR diagram is unknown. 

There are three reasons to adopt LHS 2734A as a reliable subdwarf anchor. First, we 
measure a zero parallax, implying that its Vtan is at least 280 km/sec (assuming a distance 
of 100 pc). A tangential velocity of this size is indicative of a subdwarf. Second, LHS 2734B 
{y - K, = 3.90), is redder than LHS 2734A {V - Ks = 3.23), and later type subdwarfs 
are more easily separated from dwarfs, as shown in Figure [91 As discussed in §7.2.61 LHS 
2734B is clearly a subdwarf. Assuming this common proper motion pair formed at the same 
time with similar "genetics" , both components are subdwarfs. Third, Figure [TS] shows the 
comparison between K7.0V and LHS 2734A, in which the blue end of LHS 2734A's spectrum 
is clearly "brightened". This indicates that LHS 2734A has lower metallicity than K7.0V, 
as expected for a subdwarf. 

LHS 164 This star is four full magnitudes below the main sequence line in the inset 
in Figure [TJl We find a good match at the red end between our K7.0V standard and LHS 
164's spectrum, while the blue end is elevated. We conclude that LHS 164's metallicity is 
even lower than LHS 2734A's. 

Thus, our K7.0V standard, LHS 2734A, and LHS 164 form a sequence of decreasing 
metallicity. The CaH and Ti05 band depths significantly decrease as the metallicity drops, 
matching the effects discussed in section 14.11 and Figure [3 

DEN0515-7211 and SCR 0708-4709 These two spectra are virtually identical, and 
are assigned type K7.0VI: (Figure [TJ]) because the spectra do not match either the K7.0V 
standard, LHS 164, or LHS 2734A. DEN0515— 7211 is nearly five magnitudes less luminous 
in M ks than the main sequenc e line, so it is certainly a subdwarf. SCR 0708—4709 (/x =0'.'402 



Subasavage et al.ll2005al ) has almost the same spectrum as DENOSIS- 7211 and fall s 



in the subdwarf region on the reduced proper motion diagram (see ISubasavage et al.ll2005al ). 
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1.24- MO.OVI 

We show four stars with spectral type MO.OVI in Figure [T6l This is the first group in 
this manuscript for which we see clear effects of both metallicity and gravity. 

LHS 244, LHS 418, and LHS 165 These three spectra have the same continua at 
the red end but are very different at the blue end. The differences are caused by metallicity, 
as seen in Figure [H with the main sequence standard-LHS 244-LHS 418-LHS 165 trending 
to lower metallicities. The blue ends of the subdwarf spectra have clearly "brightened" along 
this sequence, which is consistent with the theoretical models. In addition, as discussed in 
§4.11 CaHl absorption decreases (increasing index) as metallicity decreases if Tg// is greater 
than 3500K. LHS 165 is redder in V — Ks and has stronger CaHl absorption than LHS 418, 
so we suspect LHS 165 has a higher gravity (similar to gravity effects shown in Figure [T^ . 
We also find that LHS 418 has weaker CaHl absorption than LHS 244, indicating that LHS 
244 could also possibly have higher gravity than LHS 418. LHS 244 is much redder than LHS 
418 and LHS 165 in the HR diagram is because of a combination of gravity and metallicity 
effects. 

LHS 424 This star has the same spectrum and color as LHS 418, so we assign the same 
spectral type. We note that LHS 424 is a full magnitude less luminous in M^s, perhaps 
hinting that LHS 418 is a multiple system. 



LHS 300 AB This is a close binary with a separation of ~4'.'3 (jJao et al.ll2003l ) com- 



prised of a cool subdw arf (A) and a DC-type white dwarf (B) with a featureless spectrum 



(IMonteiro et al.ll2006l ). The spectrum obtain ed includes both components, but with AVRI 



= 4.61, 4.85, 4.96 mag (iMonteiro et al.ll2006l ). the contamination of the cool subdwarf spec- 
trum from the white dwarf is negligible. At a distance of 31.0 pc, the projected separation 
implies a distance of 133 AU between the two components, indicating that the subdwarf's 
composition was unlikely to be significantly contaminated by the evolved star during its 
planetary nebula phase. The spectrum does not match other MO.OVI stars at the blue end, 
but it matches the red end of MO.OV and MO.OVI well. We therefore assign it a type of 
MO.OVI:. 



7.2.5. MO. 5 VI 



This is the first half spectral type we currently assign for subdwarfs. We see both 
metallicity and gravity effects for M0.5VI, as shown in Figures [T71 and [THl 
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LHS 507, SCR 0709-4648, and LHS 12 These three stars, along with our M0.5V 
standard, illustrate a nice four-step trend in metalhcity. With the lowest metallicity, LHS 507 
is significantly bluer than LHS 12, which is, in turn, bluer than our main sequence standard. 
We do not yet have a trigonometric parallax for SCR 0709—4648, so we cannot confi rm its 



location on the HR diagram. LHS 12 was previously reported to be a type sdMO.O in iGizis 



(119971 ). 



LHS 401 This star has stronger CaH bands and therefore higher gravity than LHS 
507, which has extremely weak CaHl. LHS 401 is redder than LHS 507, as expected, but is 
slightly more luminous, which is not consistent with the gravity effects seen for types K6.0VI 
and MO.OVI. We believe this inconsistency is a result of the large absolute magnitude errors 
for these two objects. 

LHS 521 This star has a spectrum nearly identical to LHS 507, except at 6000-6200A. 
Its CaH2 feature is slightly different from LHS 507, so we assign it a type of M0.5VI:. 

LHS 367 and LHS 299 Both objects have trigonometric parallaxes. Their spectra 
are similar to SCR 0709—4648 and LHS 12, respectively (shown in Figure [TS]), so they also 
have similar offsets in metallicities from one another. We assign them as M0.5VI: because of 
(1) discrepancies at 6000A-6200A, (2) LHS 299 has deeper CaHl absorption than LHS 12 
and different CaH2 absorption, and (3) their locations on the HR diagram are redder than 
their respective comparison stars. LHS 299 may have a higher gravity than LHS 12, but a 
smaller absolute magnitude error is needed to confirm this speculation. 

LHS 360 This star appears to have higher gravity than LHS 12 because the CaH 
features are deeper, but its location on the HR diagra m does not match its stronger gravity 



because its V — Ks color is slightly bluer than LHS 12. ICosta et al.l (120061 ) reported that LHS 
360 has Vtan ~ 524 km s~^. Based on its very high Vtan and location on the HR diagram, we 
assign it a type of M0.5VI:. 



7.2.6. Ml.OVI 



The Ml.OVI type has the largest number (23) of confirmed subdwarfs in our current 
sample. The spectra of 16 different Ml.OVI stars are shown in Figures [T^ and [2U1 Many 
of these stars have similar colors, as identified in a reduced pro per motion diagram deriv ed 
using results from our SuperCOSMOS-RECONS (SCR) survey ( Subasavage et al. 2005al lbl). 
Given the rich dataset for this type, we first discuss the sample in terms of metallicity and 
gravity effects, then discuss individual targets. 
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Metallicity effects Perhaps better than for any other spectral type available, the Ml.OVI 
stars show a beautiful trend in metallicity in their spectra, as shown in the top of Figure [T9l 
The red ends of the spectra match the Ml.OV standard spectrum, but the blue ends of the 
spectra are very different because of metallicity effects. We assign their metallicities on a 

scale of m to m , where m is indistinguishable from the main sequence standard 

and 6 "-" indicate the most severely metal poor subdwarf. 

As metallicity drops, the Ti05 band gradually weakens, as predicted in the model 
spectra of Figure [HI for stars cooler than ~4000K. Our Ml.OV standard (metallicity m), 

LHS 109 (metallicity scale m ), and LHS 518 (metallicity scale m ) 

have trigonometric parallaxes, and as predicted from the models, their positions on the 
HR diagram shift bluer with decreasing metallicity. At metallicity m— , SCR 2101—5437 
is slightly metal poor compared to Ml.OV (Ti05), and has slightly higher gravity than 
Ml.OV (CaH). Unfortunately, none of the SCR objects in Figure [19] have trigonometric 
parallaxes, yet. Thus, they cannot be plotted on the HR diagram, and the progressive trend 
in metallicity effects for these stars cannot yet be shown. 

Gravity effects The effects of gravity can be seen at three different metallicities for type 

Ml.OVI. Spectral for metallicity scale m subdwarfs are shown in the middle of 

Figure [T9l SCR 0701—0655 has higher gravity than LHS 518. Both spectra are identical, 

except at CaH. Spectra for metallicity scale m subdwarfs are shown in the bottom 

of Figure dni A clear trend can be seen. Three of these five objects have trigonometric 
parallaxes — LHS 109, LHS 385 and LHS 335 — and their gravity differences shift their 
locations on the HR diagram toward redder and less luminous territory, as with previously 
discussed types. (LHS 385 's high parallax error is the likely cause of the slight offset between 
its HR diagram position and LHS 109's.) Finally, as shown at the top of Figure [20] for 
metallicity scale m , SCR 0654—7358 has higher gravity than SCR 1756—5927. 



LHS 440 was reported by iBidelmanl (Il985l ) to be type Ml.OV. The second set of spectra 



in Figure [20] shows that LHS 440 has deeper absorptions in all three CaH bands than the 
Ml.OV standard, but is virtually identical at all other wavelengths. This indicates that LHS 
440 probably has higher gravity than Ml.OV, making it less luminous than main sequence 
stars on the HR diagram. Because it is just barely one magnitude below the main sequence 
fit, we currently assign it a type of Ml.OVL. 

LHS 1970 and LHS 2734B both have noisy spectra because they are faint. As 
shown in Figure [2D1 their strong CaH bands indicate that they are high gravity (strong 
CaH), low metallicity (weak Ti05 ) subdwarfs . LHS 1970's location on the HR diagram 



confirms that it is a subdwarf and iGizisI (119971 ) reported its spectral type is esdM2.5. LHS 



2734B's spectral similarity to LHS 1970 indicates that it also has low metallicity and high 
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gravity. Unfortunately, we measure zero parallax for the LHS 2734AB pair, within the errors, 
so we can not plot them on the HR diagram (see section [7.2.31 for discussion of LHS 2734A). 
We assign both LHS 1970 and LHS 2734B types of MLOVI: because of their noisy spectra. 

LHS 158 has a spectrum nearly identical to an Ml.OV star, but falls one magnitude 



below the main sequence line on the HR diagram. It also has Vtan = 191 km sec ^ (jJao et al. 



20051 ). The low luminosity and tangential velocity together imply that LHS 158 may be a 
subdwarf, so we assign it a type of M1.0[VI], representing its uncertain assignment as a 
subdwarf, as is the case for the mid K-type subdwarfs. 



7.2.7. M2.0VI 

This is the earliest type for which we see that CaHl absorption deepens as the metallicity 
decreases, as shown by the curves in Figure [H] for T^ff less than 3500K. Figure shows the 
spectra for stars discussed in this section. 

LHS 406, GJ 191 (Kapteyn's Star) and WT0233 form a sequence of decreasing 
metallicity. These three stars have parallaxes, and their positions on the HR diagram show 
a clear trend toward the blue for stars with decreasing metallicity (shown with open arrows). 
One of the prototype s of th e subdwarf class, Kapteyn's Star, which has [Fe / H] = — . 9 9 from 



Wool 



' fc WallersteinI (120051 ) . has previously been reported to be a type sdMl.O in iGizis 



(119971 ). We assign a somewhat later type for this famous star, which at 3.9 pc is the nearest 



known subdwarf of any type. 

GJ 191/LHS 3620 and LHS 406/LHS 127 The relative effects of gravity can be 
seen in the spectra of these two pairs. LHS 3620 has much stronger gravity than GJ 191, 
while LHS 127 has slightly stronger gravity than LHS 406. On the HR diagram, LHS 3620 is 
redder and less luminous than GJ 191, as expected. LHS 162 has a spectrum (not plotted) 
virtually identical to LHS 3620, but is only slightly redder and less luminous than GJ 191. 

LHS 318 is almost identical to GJ 191 except for weaker Ti05 absorption. This weaker 
Ti05 indicates LHS 318 is slightly more metal poor than GJ 191, so the blue end of the 
spectrum should be "brightened"; however, it is not. LHS 318's location on the HR diagram 
in relation to GJ 191 suggests a higher gravity that is not apparent in its spectrum. We 
would expect LHS 318 to appear brighter and bluer than GJ 191 if metallicity is the only 
difference in these two stars, so we assign it as M2.0VL. 
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1.2.8. M3.0VI 



LHS 228/LHS 189 AB and LHS 326/LHS 398 are pairs of stars with nearly iden- 
tical spectra. As shown in Figure [221 LHS 326, LHS 228 and our M3.0V standard form a 
metallicity sequence. As in the case of M2.0V1, CaHl a bsorption deepens f or M3.0VI stars 
with decreasing metallicity. In the case of LHS189AB, I Cost a et al.l (120061 ) reported their 
separation is about 3" and their Ai? ~0.5mag. At a distance of 22.1 pc, the projected 
separation implies a distance of ~66 AU between the two components. This spectrum has 
combined both components, so we assign it type M3.0VI: until completely "clean" data can 
be acquired for the components, and their individual metallicity scales can be determined. 



WT 135 was previously identified by lHenry et al.l (120021 ) as type M2.5V. We assign its 
metallicity m— : because its Ti05 band is the same as seen in LHS 228, but the blue end 
of its spectrum is "brightened" (making it somewhat metal poor compared to LHS 228). It 
lacks a trigonometric parallax, so we cannot plot it on the HR diagram. 

LHS 272 with V -I = 4.29 is similar in color to WT 135 {V - I = 4.32) but appears 
to have slightly higher gravity (deeper CaH bands). Overall, the blue end of its spectrum 
(6000A-6300A) is comparable to WT 135, although with portions slightly brighter and 
portions slightly fainter. We therefore assign its metallicity to be the same as WT 135, at 
m— :. It is slightly redder and less luminous than LHS 228. We use both hollow (metallicity) 
and solid (gravity) arrows to indicate its location relative to LHS 228. The competition 
between metallicity and gravity effects in this case seems to indicate that gravity is the 
dominant factor because LHS 272's positio n moves to the red, rather than the blue. LHS 
272 was reported as sdM3.0 in Icizisl Jl997h . 



LHS 326, SCR 2204-3347/LHS 541 and SCR 1916-3638 form a clear sequence 
of gravity effects. Other than having a noisier spectrum (not shown), LHS 541 appears to 
be identical to SCR 2204-3347. LHS 541 's location on the HR diagram relative to LHS 326 
reflects the effects of gravity. 



7.2.9. M3.5VI 

This half type is assigned because four stars have spectra shown in Figure [23] with redder 
slopes than M3.0VI (compare to LHS 228, a M3.0VI subdwarf in the second set of spectra), 
but not as steep as M4.0VI (see next section). Conveniently, all four stars of this type have 
trigonometric parallaxes. 
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LHS 381 and LHS 144 have lower metallicities than observed for our M3.5V standard. 
Their spectra redward of 8200A match M3.5V relatively well and show the same trends 
predicted by the GAIA models (Figure [H]) for stars with temperatures of about 3200K. Their 
CaHl indices also match the trend discussed in section I4.1[ 

LHS 375 and LHS 515 have higher gravities than LHS 381 and LHS 144, respectively. 
Both obje cts ar e redd er and less luminous than their low gravity counterp arts on the HR 
diagra m. iGizid (119971 ) reported LHS 375 to be an esdM4.0 subdwarf and iReid and Gizis 
(120051 ) reported LHS 515 to be an esdM5.0 subdwarf. We classify both as M3.5VI with high 
gravities. 



7.2.10. M4.OVI to M6.0VI 

We identify only four subdwarfs with types M4.0VI to M6.0VI, only one of which, LHS 
2067A, currently has a trigonometric parallax. Because of the paucity of such objects, we 
do not yet have a sample sufficiently large to map out the effects of metallicity and gravity. 

When stars are cooler than ~3200K, their spectra redward of 7500A change radically 
between [m/H]=0.0 and —1.0 (see Figure [H]), but are rather more stable between [m/H] = 
—1.0 and —2.0. For example. Figure [2^ shows comparisons between LEHPM 3861 (M4.0VI) 
and our M3.0V and M4.0V standards. Clearly, if it s spec tral type is incorrectly assigned 



to be M3.0VI, there is excess red flux. iLodieu et al.l (120051 ) reported LEHPM 3861 to be a 
sdM6.0 subdwarf. 

LEHPM 3861 (M4.0VI), LHS 1490 (M5.0VI), LHS 334 (M6.0VI), and LHS 2067A 

(M6.0VI), whose spectra are shown in Figure [Ml are assigned subdwarf spectral types later 
than M3.5VL The types are based on fluxes redward of 7500A, with emphasis on three 
pseudo-continuum peak points at 8130, 8250 and 8840A (marked in Figure [2^ with dotted 
lines). As shown in Figure [3 for stars cooler than 3200K, the flux decreases at these three 
points as the metallicity drops from [m/H]=0.0 to —1.0 and this trend becomes even more 
prominent as T^ff decreases. Thus, these three peaks for a low metallicity subdwarf will not 
be brighter than a dwarf with a comparable spectral type. 

LHS 2067A is bluer than M6.0V on the HR diagram and is slightl y below the main 



sequen ce line. LHS 2067A was previously identified to be a subdwarf in iKirkpatrick et al. 



(119951 ). but no spectral type was given. Although its CaHl and CaH2+CaH3 indices are 
not in the subdwarf region (see Figure [3] and Table [1]), its spectrum is clearly different from 
M6.0V and M7.0V, with stronger CaHl absorption. It forms a wide (~55" NE) common 
proper motion pair with a white dwarf, LHS 2067B, so it, like LHS 193AB and LHS 300AB, 
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forms an intriguing pair that can be used for comparing metallicity and white dwarf ages. 
At the system's distance of 25.5 pc, the large separation of the pair corresponds to 472 AU, 
indicating that significant pollution of the subdwarf by the white dwarf during its planetary 
nebula phase seems unlikely. 



LHS 334 was reported to be a sdM4.5 subdwarf in iReid and Gizid (120051 ). We assign 
it a type of M6.0VI because of the good match of the three pseudo-continuum peak points 
to LHS 2067A, which appears to be slightly more metal rich than LHS 334. 



8. Application to SDSS Subdwarfs 



It is useful to apply our spectral typing methodology to the recent work of IWest et al. 



( I2OO4J ). who have provided a significant sample of sixty new subdwarf spectra, all acquired 
and reduced in a homogeneous way. The subdwarfs were selected and spectral types were 
assigned via their CaH and TiO indices. For comparison, we retrieved fift}@ of the publicly 
available spectra and use our new method to assign types. The brightest star among these 
subdwarfs has r=17.0, resulting in somewhat noisy spectra for the sample, so we smoothed 
the spectra by averaging the flux over five pixels and normalizing at 7500A. Representative 
spectra are shown in Figures [SS] and [ 



The SDSS s p ectra have telluric lines removed, so we use M dwarf standard spectra from 
Bochanski et al.l (120071) th a t also omit the telluric lines, rather than our own standards. In 
addition, iBochanski et al.l (120071 ) do not present half-type standard spectra, so we do not 
present any half-type spectra for the SDSS subdwarfs. Previous an d current types ar e liste d 
in Table [31 The subdwarfs have types between Ml.O and M3.0 in iBochanski et al.l (120071 ). 
while our efforts yield types between M2.0 and M5.0. 

Our point here is not to assign definitive spectral types, but to illustrate the gravity 
and metallicity effects in the SDSS dataset, and provide examples of the application of our 
proposed spectral typing method. Stars of a given spectral type, of course, may have a range 
of metalhcities, as discussed at length in previous sections, and as is evident in Figure [251 
Gravity effects are also seen, as shown in Figure [261 Three stars (SDSSJ085843. 89+511210.1, 
SDSSJ093141. 85+453914.5 and 

SDSSJ145447.32+011006.8) previously identified as subdwarfs have spectra identical to SDSS 
M dwarfs. Two other stars, SDSSJ083217.77+522408.2 and SDSSJ113501. 76+033720.3, ap- 



^The remaining 10 subdwarfs could not be retrieved using the SDSS DR4 website, even though we used 
a 1' search radius and coordinates from their table. 
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pear to have slightly higher gravities (stronger CaHl line) than dwarfs, but the rest of these 
two spectra are almost identical to SDSS M dwarfs. Therefore, to these two stars we assign 
types of M2.0VI: and M3.0VI: with g+:, respectively. 

Unfortunately, the SDSS dwarfs do not have parallaxes that can be used to confirm 
their locations on the HR diagram. In addition, we find many SDSS subdwarfs with identical 
spectral types that have g — z colors as different as 0.53. It is therefore difficult to verify how 
metallicities and gravities af f ect th e SDSS absolute magnitudes and colors. Nonetheless, the 
robust sample of I West et al.l (120041 ) provides many additional subdwarfs that can be targeted 
for further work. 



9. Discussion 

9.1. Why Subdwarfs' Physical Parameters Are Not Listed Explicitly 

Stellar spectra follow trends primarily defined by temperature, yet as we have seen, 
metallicities and gravities also have significant effects on the spectra of cool subdwarfs. 
Here we have provided a consistent spectral sequence for subdwarfs that is based on linking 
observed subdwarf spectra to spectra acquired for main sequence dwarfs using the same 
telescope/instrument /observing protocols. Here we compare synthetic to observed spectra 
to evaluate how well we can assign values for the temperatures, metallicities, and gravities 
of cool dwarfs and subdwarfs. 



9.1.1. Matching Synthetic and Observed Spectra 

To test the reliability of the model grids and our fitting procedures, we use our six 
standard spectral sequence dwarfs with types MO.OV to M5.0V as test spectra (see Fig- 
ure [27]). All are currently believed to be uncorrupted single red dwarfs because they show 
no evidence of multiplicity from combinations of (1) HST/NICMOS observations, (2) optical 
speckle observations, (3) optical CCD imaging, and/or (4) three or more years of astromet- 
ric observations that would reveal perturbations from unseen companions that contributed 
significant light to the system. 

We compared the spectral region from 6OOOA to 9000A after normalizing both the 
GAIA grid spectra and ours at 7500A, and applying a Gaussian function to the much higher 
resolution GAIA model spectra to match our resolution. If a star had a significant spectral 
shift because of radial velocity, we manually offset the spectrum to match features obvious 
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in the synthetic spectra. We then calculated the reduced differences (hereafter, simply 
X^) between the model grid spectra and ours. Because the synthetic spectra do not have 
the telluric lines of Oaa (6270A-6330A), O2B (6860A-6980A), O2A (7590A-7710A) and 
water (7150A-7330A and 8952A-9000A), these absorption regions were excluded in the 
calculations. 

The selected model grids have effective temperatures of 2400K to 4500K, [m/H] from 
—2.0 to +0.5, and log g from 4.0 to 5.5 in steps of lOOK, 0.5 dex, and 0.5 dex, respectively. 
The upper panel of Figure [TD] shows that for a metallicity of —1.0 and temperature 3500K, 
varying the gravity causes changes in only certain wavelength regions (and some specific 
lines). This means is only sensitive to log g in relatively small spectral regions, but not 
for the overall spectrum. On the other hand, as shown in the lower panel of Figure [101 for 
the same temperature star with log g fixed at 5.0, there are large differences between spectra 
when the metallicity is varied. 

Thus, for a star of a given temperature, the gravity changes the overall spectral shape 
minimally while the metallicity changes it a great deal, so we first secure a star's metallicity 
and then its gravity. Figure [251 is an example for Ml.OV standard star, showing curves 
at various metallicities and gravities. Each curve represents a specific log g and grid spectra 
have temperatures incremented by lOOK. The smallest scatter in the plots is found when 
[m/H] is —0.5, which is adopted as the star's metallicity. We then examine each point in the 
[m/H] =—0.5 panel to find the best fit with the smallest this fixed metallicity. In this 

case, the best fit has T^jf = 3600K, log g = 4.5 and [m/H] = —0.5. 

9.1.2. Discrepancies Between the Best Fitting Synthetic and Observed Spectra 

In this paper we do not explicitly list temperatures, metallicities, or log g values for 
subdwarfs. A few examples support our reasoning for not listing these physical parameters. 
Figure shows observed spectra and the best fit synthetic spectra for six main sequence 
spectral standard stars with types MO.OV to M5.0V (in steps of 1.0 subtypes). In general, the 
overall slopes of the model spectra fit fairly well, especially for the earlier types. However, 
we discuss here six regions labeled at the top of Figure [27] that do not match, which are 
particularly relevant for the subdwarfs that are the focus of this paper. 

1. Region 1: The observed spectrum is always less luminous than the model, except for 
MO.OV. 

2. Region 2: This is the CaHl absorption region. Observed spectra have shallower CaHl 
features than the models. The depth of CaHl is determined by a combination of gravity 
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and metallicity. Decreasing CaHl absorption can be caused by either decreasing gravity 
or decreasing metallicity. However, fine tuning the gravity or metallicity affects not 
only the CaHl feature, but the Ti05 feature in region 4. 

3. Region 3: The pseudo-continuum always peaks near 6530 A in the observed spectra. 
However, the synthetic spectra have this peak "red-shifted" to ~6650A. 

4. Region 4: Containing the CaH2, CaH3 and Ti05 features, this is the most important 
region for examining the interplay of metallicity and gravity in cool dwarfs. Unfortu- 
nately, the CaH2 feature is blended with O2B so it is not entirely reliable for analysis of 
spectra taken through the Earth's atmosphere. The CaH3 and Ti05 features are usu- 
ally weaker in observed spectra than in the models. The strength of Ti05 is primarily 
driven by metallicity, not gravity (see the bottom panel of Figure [TOj) . Therefore, this 
region reveals valuable information about a star's metallicity, in particular at types 
later than M2.0. 

5. Region 5: Several Fe I absorption lines (8388A, 8440A and 8718A) and a Mg I line 
(8718A) in the model spectra redward of 8300A are deeper than observed. This indi- 
cates poor metallicity matches and/or poor modeling of those particular lines. 

6. Region 6: Overall, redward of 8000A the M4.0V and M5.0V matches are not as good as 
other regions. This likely indicates fundamental problems with the strengths of some 
opacity sources (e.g. H2O and TiO) in the models. 



The top panel of Figure [291 shows the two best fitting synthetic spectra for our Ml.OV 
standard. The red spectrum (3600K, [m/H]= —0.5, log g=A.5) provides the best fit and 
the blue spectrum (3400K, [m/H]= —1.0, log g=A.O) is second best. There are only slight 
differences between the two synthetic spectra — the CaHl strengths and continuum fluxes 
redward of 8OOOA are the only notable differences. However, the effective temperatures for 
the two models differ by 200K, while the metallicities and log g each differ by 0.5 dex. The 
lower panel of Figure [22] shows the two best fit spectra for our M3.0V standard. The best 
fit (red line) yields 3200K, [m/H]= —0.5, log g=4:.5 and the second best fit (blue line) yields 
3300K, [m/H]=0.5, log g=5.5. However, neither model fit is an ideal match to the observed 
spectrum (which is why the [m/H] value for the two best matches differs by an order of 
magnitude), and matches become even poorer for cooler stars. 

As a whole, many of the best fits for the main sequence dwarfs in Figure [27] are for 
metallicities of —0.5, which is somewhat lower than studies that have specifically attempted 
to assign metallicities to cool dwarfs i n the solar n e ighbo rhood. For example, the mean 



[Fe/H] from 21 M dwarf secondaries in [Bonfils et al.l (120051 ) is —0.09, and the mean [Fe/H 
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from five M dwarf secondaries in lBean et al.l (120061 ) is —0.17. Also worthy of note is that the 



fits were made for metallicities incremented by 0.5 dex, and stars with metalhcities between 
—0.5 and 0.0 may have slightly better fits for —0.5. One might think that interpolation 
from existing model grids would allow better fits. There are, however, degeneracy problems 
in matching model spectra to observations — when fitting the three-dimensional space of 
temperature, metallicity, and gravity, two or more synthetic spectra yield low points with 
similar , as is the case for the Ml.OV example discussed above. More worrisome is that the 
overall discrepancies discussed above (items 1 to 6) cannot be removed through interpolation 
in the existing model grids. 

Figure [30] illustrates how the values change for best fit matches of the GAIA models to 
our spectral standards with types MO.OV through M5.5V. Spectral type M3.0V is the latest 
type for which a reasonably tight plot like the one shown in Figure [28] can be identified. 
The dotted line in Figure [30] provides a dividing point between reasonable matches and poor 
matches, indicating that once exceeds 10, the plots are too scattered to choose a unique 
set of model spectra parameters to match observed spectra. 

A specific example of applying the models to one of our subdwarfs illustrates the large 
discrepancies that must be overcome to derive reliable parameter values for metallicity and 
gravity. When we applied our fitting algorithm to our observed spectrum for LHS 335 
(Ml.OVI, shown in Figure [31]), we discovered that the "best fit" was quite poor. Although 
the formal value was 5.9 (less than our cutoff of 10, indicating a reliable fit), the model's 
CaH3 band is not deep enough, the Ti05 band is too deep, and the continuum flux redward 
of 7700A is less than observed. 

For myriad reasons we conclude that we cannot strictly determine reliable metallicities 
and log g values for cool dwarfs using the fitting method discussed here. Thus, until im- 
proved model grids are available, we defer assignment of numerical values for temperatures, 
metallicities and log g values. However, we certainly can use the GAIA model grids to mimic 
trends (discussed in §5] and [6]) in the spectra of subdwarfs to compare the stars within a 
framework of changing physical parameters. 



9.2. The Confusion Between sd, esd and usd Prefixes 



Gizid (119971 ) proposed that subdwarfs have two subclasses, "subdwarfs" (sd) and "ex- 
treme subdwarfs" (esd) based on their CaH and Ti05 band strengths. He separated the two 
classes using the line shown in the top panel of Figure [3] Another term, "ultra subdwarf" , 
was proposed bv ICaldwell et al.l ^M) to describe GJ 59B {V - K, = 3.42, Mks = 7.89), 
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which is ~2.5 magnitudes underluminous compared to main se quenc e stars of similar color. 
This star fall s in th e region including "extreme subdwarfs" in iGizij (Il997l ). More recently, 



Lepine et al.l (120071 ) adopted the "ultra subdwarf" (usd) term for subdwarfs found to have 
stronger CaH features than "extreme subdwarfs" . We believe these terms confuse the situa- 
tion and do not address the underlying astrophysics. We recommend that they not be used 
for the following reasons. 



• These terms can only be applied to M-type subdwarfs, not K-type subdwarfs. Both 
observed spectra and models show that K subdwarfs do not follow the trends in CaH 
absorption with metallicity that M subdwarfs follow. Thus, there is no clear delineation 
for K-type subdwarfs. 

• Empirically, the values of the indices are affected by a complicated interplay of temper- 
ature, metallicity, and gravity effects. One cannot separate these three factors simply 
based on the indices. Typically, when the term "extreme" subdwarf is used, it refers to 
"very low metallicity" alone. But that is only one part of the portrait that needs to be 
painted for a given subdwarf. For example, filled boxes shown in Figure H] indicate stars 
previously identified as extreme subdwarfs from their indices. Note that the spread for 
"extreme" subdwarfs at a given color (V — Kg 3.5) is 3.5 magnitudes, or a factor 
of ~25 in luminosity. There is no clear separation in the fundamental HR diagram 
between stars termed "subdwarfs" and "extreme subdwarfs" classified using spectral 
indices. 

In Figure [321 we take a detailed look at astrophysical causes that shift points in the 
CaH2-|-CaH3/Ti05 indices plot. Red circles in Figure [32] represent the lowest metal- 
licity stars of each spectral type we have presented here. Some of these very low 
metallicity stars are not located in the "extreme" subdwarf region. Green circles rep- 
resent subdwarfs in our sample that have highest gravity at a given type. Solid lines 
connect stars having the same metallicity rankings in Table [2] (components of the com- 
mon proper motion binary LHS 2734 AB are connected by a dotted line because they 
presumably have the same metallicity). It is clear from these three pairs that high 
gravity can push a subdwarf toward or into the extreme subdwarf region, even if the 
metallicities of the two objects are similar. In addition, many other high gravity sub- 
dwarfs (green circles) that do not have particularly low metallicities are also in the 
subdwarf region. Finally, the two blue circles represent LHS 440 and SCR 1822-45542, 
which subdwarfs having the same metallicities but higher gravities than dwarfs. Their 
spectroscopic features indicate that they are subdwarfs, but they do not have low 
metallicity at all. Yet, they fall in the subdwarf region of the indices plot, which has 
traditionally indicated low metallicity. 
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The inset in Figure [32] helps explain this phenomenon. We have calculated the spec- 
tral indices from the GAIA synthetic spectra with [m/H]=0.0, —1.0 and —2.0 and 
temperatures of 2800-4400K. Two different gravities were selected for each metallicity. 
The models indicate that both low metallicity and high gravity push stars toward the 
extreme subdwarf region. Thus, if two stars have the same metallicity, the one with 
higher gravity will be pushed more toward the extreme subdwarf region, indicating 
that this spectral indices plot is not a clear indicator of metallicity alone. 

• Theoretically, in order to show how complicated the Ti05 vs CaHl plot is astrophysi- 
cally, we use GAIA model grids to do a demonstration. We have chosen models with 
2700K < Teff < 4500K, -3.0 < [m/H] < 0.0 and 4.0 < logg < 5.5, and have calculated 
the output CaH and Ti05 indices, as shown in Figure [331 As is apparent in the plots, 
for a given (Ti05, CaHl) indices pair, there are many possible parameter combina- 
tions. Both indices are a function of temperature, gravity and metallicity, which would 
require a 6-D plot to describe Ti05 vs CaH. Unfortunately, one cannot classify stars 
as subdwarfs, extreme subdwarfs, or ultra-subdwarfs using 2-D index plots that infer 
trends in metallicity alone. 

For these reasons, the indices only indicate a star's location on these particular spec- 
troscopic indices plots. They do not address low or very low metallicity only, as implied 
by past usage, but instead incorporate effects of both metallicity and gravity for stars of 
a given temperature. They do not provide obvious or direct links to positions on HR di- 
agrams, which provide the astrophysical meaning underlying stellar classification. Finally, 
the modifiers "extreme" and "ultra" themselves have effectively the same meaning and do 
not provide meaningful information about their differences. 



9.3. Why Previous Methods Work for Dwarfs but Not Subdwarfs 



In the Palomar-MSU spectroscopy survey, iReid et al.l ( 19951) us ed the Ti05 index and 
polynomial equations to assign subtypes for M dwarfs. iGizisI (119971 ) then applied the same 
methodology to M subdwarfs. As can be seen in Figure [331 Ti05 versus temperature is 
effectively linear, especially for log g = 5.0 and [m/H] = 0.0, so this typing method works 
successfully for a fairly homogeneous set of dwarfs (possibly with somewhat different grav- 
ities) — temperature is the main factor affecting the overall slope of red dwarf spectra. 
However, it is not as simple for subdwarfs because both low metallicity and high grav- 
ity stars are sorted into the same Ti05 values for various combinations. Thus, the same 
methodology that works for dwarfs cannot straightforwardly be applied to subdwarfs. 
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9.4. The Subdwarf Spectral Standards from iLepine et al.l (120071 ) 



Lepine et al.l (120071 ) recently released a set of spectral stan dards f or the sdM, esdM 



and usdM subclasses using wider spectral coverage than used by iGizij ( 119971 ) and revised 
polynomial equations. Their project is contemporaneous with ours, but unfortunately there 
is only one star in both samples, LHS 228. We assign it a type of M3.0VI while their type is 
sdM2.0. As shown in Figure [311 gravity effects come into play in their spectre^. Each pair 
of spectra shown have nearly identical continua, with CaH bands being the only significant 
spectral difference. Using their method, each of the pairs of subdwarfs shown is assigned 
a different sub-type. We believe that by using the overall shape of the spectrum, with 
knowledge that both metallicity and gravity affect certain regions, each pair should have the 
same spectral sub-types. Based on the trends from GAIA models, the pairs of spectra shown 
have different gravities, as we have also seen in our spectra and those from SDSS (shown in 
Figure [35] illustrates spectral differences at types M3.0, M5.0 and M7.5. Although the 
spectra in each panel have the same sub-types for sdM, esdM and usdM subdwarfs, they 
show very different continua between 6OOOA and 9000A. When examining the overall spectra, 
there is no clear morphology trend for each sub-type. We conclude that using spectral indices 
alone omits important information evident in the overall morphology of the spectra that is 
useful to spectral classification. 



10. Conclusions 

We have discussed 88 cool subdwarfs using spectra covering 6OOO-9OOOA. Based on these 
spectra and the trends from GAIA model grids, we have redefined the subdwarf spectral 
sequence, spanning types K3.0 to M6.0. We find that wide spectral coverage is the key to 
defining a subdwarfs spectral type. We consider this to be an important, but not final, step 
in defining the subdwarf spectral sequence. 

Through the understanding of GAIA model grids, we find that the key to assigning a 
subdwarf spectral type is to compare the spectrum to dwarf spectral standards in regions 
affected minimally by metallicity and gravity, thereby making a direct link between the dwarf 
and subdwarf sequences. Even so, it remains difficult to establish a definitive sequence for 
subdwarfs because of the multi-faceted nature of their spectra. Until we have surveyed a 
large number of subdwarfs and covered a multitude of possible temperatures, metallicities, 
and gravities, a definitive sequence will remain elusive. 



Their telluric lines have been removed. 
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From an analysis of the history of the term "subdwarf," and the layout of the funda- 
mental HR diagram, we propose that the suffix "VI" be used, rather than the "sd" prefix, 
as the preferred spectral classification notation. This reduces the confusion between cool 
subdwarfs and hot OB subdwarfs. This is also prudent because subdwarfs really do form 
an independent class of stars on the HR diagram, for which five Roman numerals are in 
common use, with cool subdwarfs naturally falling beneath the main sequence V types. 

Overall, we find that trigonometric parallaxes are crucial for identifying mid K-type 
subdwarfs, and allow us to understand how the complex interplay of temperatures, metal- 
licities, and gravities affects the positions of cool dwarfs on the HR diagram. We found 
that mid K-type subdwarfs can not be identified spectroscopically using our data. There 
are many G and M-type subdwarfs, so K-type subdwarfs presumably exist, but additional 
information is needed to identify them. Consequently, we use metallicities, kinematics, and 
parallaxes to acertain their true natures and use their locations on the HR diagram to fiag 
them as "subdwarfs" (other spectroscopic wavelength and resolution combinations could also 
be used). Because they are not spectroscopically identified as subdwarfs in our spectra, we 
use a conservative notation, [VI], for subdwarfs of types K3.0 through K5.0 to indicate their 
questionable status. 

We have confirmed that spectroscopic indices are useful in separating late K to late 
M-type subdwarfs from dwarfs, but that the indices have limitations when attempting to 
understand the astrophysical causes leading to observed subdwarf spectra. When combined 
with trigonometric parallax and photometric information, our results show that for lower 
metallicities, subdwarfs are generally bluer and brighter at optical wavelengths, so they 
slide up and to the left on the HR diagram (using axes V — Ks vs Mks)- In contrast, 
higher gravities make stars redder and less luminous at optical wavelengths, so subdwarfs 
generally slide down and to the right on the same HR diagram. Because of the complex, 
and not yet completely mapped out, interplay of temperatures, metallicities, and gravities, 
we conclude that the "extreme" and "ultra" prefixes only outline locations on spectroscopic 
indices plots, and do not successfully differentiate the underlying astrophysical causes for 
shifts on CaH/Ti05 plots. 

Improvements in the subdwarf spectral sequence can be made by observing wide common 
proper motion subdwarf binaries, like LHS 2734AB discussed here. Assuming identical 
metallicities, such binaries allow us to constrain one of the three variables that affects cool 
dwarf spectral types. The ultimate subdwarf spectral sequence will be three-dimensional 



with temperature, metallicity, and gravity as independent variables (see iKirkpatrickl 12005 
Figure 11). 

Previously, gravity effects in cool subdwarfs have been almost entirely ignored, as metal- 



- 35 - 



licity was the factor considered to describe changing CaH features. Our results clearly show 
the importance of gravity effects. An additional factor, contamination from an unseen, 
evolved, i.e. white dwarf, companion could also change the slope of subdwarf spectra. Thus, 
comprehensive surveys of subdwarfs for companions are warrantee!^. Finally, this work also 
shows that current synthetic spectra provide a useful framework in which to evaluate cool 
dwarf spectra, but they do not yet provide perfect matches, indicating that atmospheric 
models still require fine-tuning to make additional advances in the future characterization of 
cool dwarfs (P. Hauschildt private communication). 
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Table 1. Spectroscopic Indices 
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.033 


LHS2734A 




1. 


,018 


0. 


,985 


0. 


,984 


0, 


.987 


1 


.972 


LHS2734B 




1. 


,086 


0. 


,720 


0. 


,669 


0, 


.828 


1 


.497 


LHS3620 




0. 


,836 


0. 


,684 


0. 


,544 


0, 


.727 


1 


.271 


SCR0242- 


5935 


0. 


,910 


0. 


,906 


0. 


,838 


0, 


.891 


1 


.728 


SCR0406- 


6735 


1. 


,015 


0. 


,989 


0. 


,829 


0, 


.937 


1 


.766 


SCR0433- 


7740 


0. 


,929 


0. 


,887 


0. 


,760 


0, 


.864 


1 


.624 


SCR0529- 


3950 


0. 


,817 


0. 


,886 


0. 


,637 


0, 


.815 


1 


.452 
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Table 1 — Continued 



01)ject 


Ti05 


CaHl 


CaH2 


CaH3 


CaH2 


+CaH3 


SCR0629- 


-6938 


0. 


o A 

,845 


0. 


,908 


0. 


,579 


0, 


,784 


1, 


.363 


bCR0654- 


-7358 


0. 


,938 


0. 


o o r 

,835 


0. 


,653 


0, 


.825 


1 


.478 


SCROyOI- 


-0655 


0. 


,976 


0. 


,886 


0. 


,780 


0, 


.881 


1 


.662 


bCR0708- 


-4709 


0. 


,984 


0. 


,992 


0. 


,979 


0, 


.976 


1 


.955 


SCR0709- 


A A a 

-4648 


1. 


,013 


0. 


,899 


0. 


,846 


0, 


.915 


1 


.761 


bCR1107- 


-4135 


0. 


,994 


0. 


o o r\ 

,839 


0. 


,783 


0, 


.892 


1 


.675 


SCR1433- 


-3847 


1. 


,001 


0. 


,904 


0. 


,797 


0, 


.930 


1 


.726 


SCR1455- 


-3914 


0. 


,915 


0. 


'"70 o 

,788 


0, 


,707 


0, 


.849 


1 


.555 


SCR1457- 


-3904 


0. 


,905 


0. 


,824 


0. 


'^r\ '~7 

,707 


0, 


o o o 

.838 


1 


.545 


bCR1613- 


-3040 


0. 


,920 


0. 


,898 


0. 


,749 


0, 


o o 

.885 


1, 


.634 


bCR1739- 


-8222 


0. 


,945 


1. 


,009 


0. 


,797 


0, 


.882 


1 


.679 


bCR1740- 


-5646 


0. 


,915 


0. 


,696 


0. 


,496 


0, 


.637 


1 


.133 


SCR1756- 


-5927 


0. 


,922 


0. 


,794 


0. 


,757 


0, 


.856 


1 


.613 


SCR1822- 


-4542 


0. 


,817 


0. 


,818 


0. 


,624 


0, 


.809 


1 


.433 


SCR1843- 


-7849 


0. 


,869 


0. 


,810 


0. 


,715 


0, 


.829 


1 




SCR1913 


-1001 


0. 


,962 


0. 


,923 


0. 


,835 


0, 


.925 


1 


.760 


SCR1916- 


-3638 


0. 


,931 


0. 


,647 


0. 


,493 


0, 


.666 


1 


.159 


SCR1958- 


-5609 


0. 


,915 


0. 


,865 


0. 


,805 


0, 


.926 


1 


.731 


SCR2018- 


-6606 


0. 


,886 


0. 


,856 


0. 


,715 


0, 


.836 


1 


.550 


SCR2101- 


-5437 


0. 


,929 


0. 


,969 


0. 


,726 


0, 


.904 


1 


.630 


SCR2104- 


-5229 


0. 


,940 


0. 


,907 


0. 


,772 


0, 


.865 


1 


.638 


SCR2109- 


-5226 


0. 


,967 


0. 


,774 


0. 


,646 


0, 


.799 


1 


.445 


SCR2204- 


-3347 


0. 


,874 


0. 


,711 


0. 


,569 


0, 


.762 


1 


.331 


SIP1342- 


3534 


0, 


,684 


0, 


,664 


0, 


,462 


0, 


.676 


1 


.138 


WT0135 




0, 


,622 


0, 


,754 


0, 


,512 


0, 


.731 


1 


.243 


WT0233 




0. 


,939 


0. 


,933 


0. 


,739 


0, 


.857 


1 


.596 



Table 2. Spectral lypes 



RA 


DEC 


Object 




V - Ks 


^ s 


Vtan 


Old Type 


New Type 


metallicity 


gravity 




Rcf 
















km /sec 










7r 


V 


spcct 


12 06 00.9 


+14 38 56.8 


G012-016 


7.931+0.018 


2.16 


4.12 


98.4 


K2 


GVI 


(-0.52)^ 




H 


2 


2 


01 04 26.4 


-02 21 59.8 


G070-035 


7.116+0.020 


2.04 


4.15 


55.8 


G5 


GVI 


(-0.67)^ 




H 


4 


H 


02 25 49.8 


+05 53 39.5 


G073-056 


10.427+0.023 


2.06 


8.64 


44.9 




GVI 


(-1.18)'' 




H 


4 




U7 54 ii.l 


— Ul 24: 44. o 


Cjrll2— 054 


5.425+0.023 


2.00 


4.01 


23.7 


Kl 


GVi 


(—0.94)^ 




H 


1 


2 


lb 1.^ 48. o 


— 57 34 13.8 


T TUOn /I 1 o 

Lllb0413 


5.293+0.024 


2.24 


4.60 


106.5 


G8/K0V(Wj 


GVI 


(-1.35)'= 




H 


1 


10 


15 45 52.4 


+05 02 26.6 


G016-009AB 


6.880+0.024 


2.27 


3.46 


59.2 


K2V 


K3.0[VI] 


m-, (-0.77)=^ 


g 


H 


H 


3 


00 50 17.0 


-39 30 08.3 


LHS0125 


11.452+0.026 


2.88 


7.73 


342.6 




K4.0[VI] 


m— 


g 


R 


R 




07 35 46.3 


+03 29 36.0 


LHS0232 


10.841+0.024 


2.84 


6.57 


346.0 




K4.0[VI] 


m— 


g 


Y 


R 




12 25 50.7 


-24 33 17.8 


LHS0327 


10.144+0.021 


2.57 


5.44 


418.3 


KOV 


K4.0[VI] 


m— 


g 


R 


R 


2 


11 52 32.0 


+27 .30 51.3 


LHS2467 


9.806+0.017 


2.43 


4.94 


437.0 


G7V 


K4.0[VI] 


III— 


g 


H 


R 


2 


19 07 02.0 


+07 36 57.3 


G022-015 


6.469+0.018 


2.71 


4.54 


97.2 


K5V 


K5.0[VI] 


m-, (-0.61)'* 


g 


H 


H 


2 


21 32 11.9 


+00 13 18.0 


G026-009ACD 


7.082+0.029 


2.64 


3.62 


96.9 


K2V 


K5.0[VI] 


m-, (-1.05)f 


g 


H 


H 


2 


05 54 34.1 


-09 23 33.7 


GJ0223.1 


7.760+0.020 


2.96 


4.84 


82.1 


K4V 


K5.0[VI] 


m-, (-0.62)"= 


g 


H 


1 


9 1 


04 32 36.6 


-39 02 03.4 


LHS0193A 


8.427+0.023 


3.23 


5.80 


147.3 




K6.0VI 


m— 


g 


R 


16 




23 43 16.7 


-24 11 16.4 


LHS0073 


9.393+0.021 


3.42 


7.27 


322.3 


K5V 


K6.0VI 


m— 


g+ 


Y 


19 


2| 


07 13 40.6 


-13 27 57.1 


LHS0227 


11.036+0.023 


3.41 


7.15 


362.5 




K6.0VI 


m— 


g+ 


Y 


R 




02 52 45.7 


+01 55 49.4 


LHS0161 


10.995+0.019 


3.65 


8.01 


272.6 


esdM2 


K6.0VI 


m— 


g++ 


Y 


R 


8 


03 01 40.6 


-34 57 56.5 


LHS0164 


10.641+0.024 


2.92 


8.20 


192.6 




K7.0VI 


m 


g 


R 


R 




05 15 45.1 


-72 11 22.2 


DEN0515-7211 


13.211+0.040 


3.29 


9.72 


4.7 




K7.0VI: 


m 


g 


7 


7 




07 08 32.0 


-47 09 30.5 


SCR0708-4709 


10.764+0.025 


3.05 








K7.0VI: 


m 


g 




R 




13 25 14.0 


-21 27 06.0 


LHS2734A 


12.896+0.037 


3.23 


too far 


>280 




K7.0VI 


m— 


g 


R 


R 




03 06 28.7 


-07 40 41.5 


LHS0165 


11.006+0.025 


3.42 


7.67 


333.1 




MO.OVI 


m 


g+: 


Y 


R 




11 11 13.7 


-41 05 32.7 


LHS0300A 


9.802+0.023 


3.38 


7.23 


197.5 




MO.OVI: 


m 


g 


11 


16 




16 25 14.0 


+15 40 54.2 


LHS0418 


10.072+0.018 


3.37 


6.44 


301.5 


K7V 


MO.OVI 


m 


g 


Y 


R 


8 


16 37 05.6 


-01 32 01.6 


LHS0424 


10.803+0.022 


3.37 


7.39 


279.6 




MO.OVI 


m 


g 


Y 


R 




08 13 27.8 


-09 27 56.6 


LHS0244 


10.729+0.023 


3.64 


6.88 


412.7 




MO.OVI 


m— 


g++: 


Y 


19 




22 27 59.0 


-30 09 30.0 


LHS0521 


11.463+0.019 


3.22 


8.14 


222.1 




M0.5VI: 


m 


g 


11 


11 




21 21 34.8 


-19 03 38.6 


LHS0507 


12.053+0.026 


3.15 


7.81 


355.2 


K/M sd 


M0.5VI 


m 


g+ 


Y 


R 


Y 


15 39 39.0 


-55 09 10.0 


LHS0401 


9.407+0.019 


3.31 


7.33 


142.2 




M0.5VI 


m 


g++ 


Y 


R 




02 42 26.3 


-59 35 01.6 


SCR0242-5935 


12.783+0.031 


3.40 








M0.5VI 


m 


g 




R 




07 09 37.2 


-46 48 58.8 


SCR0709-4648 


11.491+0.026 


3.41 








M0.5VI 


m 


g 




R 




11 07 55.8 


-41 35 52.7 


SCR1107-4135 


11.474+0.019 


3.49 








M0.5VI 


m 


g 




R 




17 39 45.4 


-82 22 02.2 


SCR1739-8222 


12.190+0.026 


3.48 








M0.5VI 


m 


g 




R, 




02 02 52.2 


+05 42 21.0 


LHS0012 


8.684+0.020 


3.56 


6.32 


342.1 


sdMO 


M0.5VI 


m— 


g 


Y 


13 


8 



Table 2 — Continued 



RA DEC Object Ks V - Ks Mk, Vtan Old Type New Type metallicity gravity Ref 

km/sec tt V spect 



11 11 22.6 -06 31 56.4 LHS0299 11.143±0.024 3.64 6.54 437.3 M0.5VI: m- g Y R 
14 18 20.4 -52 24 12.6 LHS0367 9.786±0.019 3.11 6.26 268.6 M0.5VI: m- g 7 7 
04 06 06.7 -67 35 28.8 SCR0406-6735 12.804±0.031 3.54 •■• ••■ M0.5VI m- g R 

04 33 26.5 -77 40 09.7 SCR0433-7740 13.361±0.034 3.50 •■• ••• M0.5VI m- g R 
14 33 03.3 -38 46 59.6 SCR1433-3847 13.592±0.046 3.62 • ■ • • • ■ M0.5VI m- g R 

18 43 35.7 -78 49 02.5 SCR1843-7849 12.591±0.026 3.55 • • • • • • M0.5VI m- g R 

19 58 31.2 -56 09 10.6 SCR1958-5609 12.525±0.033 3.51 • ■ • • • ■ M0.5VI m- g R 

13 46 55.5 +05 42 56.4 LHS0360 11.662ili0.023 3.50 6.73 524.0 M0.5VI: m- g-j- 7 7 

22 20 27.0 -24 21 49.3 LHS0518 10.393±0.021 3^24 TOT 230.7 Ml.OVI m i 11 U 

07 01 17.7 -06 55 49.3 SCR0701-0655 12.996ili0.030 — — Ml.OVI m g+ 

01 53 09.0 -33 25 02.1 LHS0148 12.832±0.032 3.59 8.59 374.8 Ml.OVI m g+: 6 6 

01 53 09.0 -33 25 02.1 SCR1913-1001 11.929±0.028 3^69 ~ ~ Ml.OVI m i R 

00 17 40.0 -10 46 16.9 LHS0109 10.366±0.021 3.51 7.63 176.1 K5 Ml.OVI m g+ Y R 2 

21 04 00.5 -52 29 43.5 SCR2104-5229 12.763±0.024 3.53 ••■ ■■■ Ml.OVI m g+ R ' 

14 55 35.8 -15 33 44.0 LHS0385 11.062±0.023 3.55 7.61 403.4 MO Ml.OVI m g++ Y R Y^ 

14 55 51.5 -39 14 33.1 SCR1455 -3914 11.788±0.024 3.65 ••• ••• Ml.OVI m g++ R | 

12 34 53.1 +05 03 54.1 LHS0335 13.000+0.030 3.60 9.74 247.2 Ml.OVI m g+++ Y R 

21 09 02.5 -52 26 17.8 SCR2109-5226 13.049+0.036 3.66 • • • • • • Ml.OVI m g+++ R 

01 31 04.1 -50 24 54.3 LEHPM1628 13.459+0.038 3.70 •■■ ■■• Ml.OVI m g++++ R 

06 54 06.3 -73 58 03.6 SCR0654-7358 13.285+0.048 ~ ~ ~ Ml.OVI m i 

14 57 49.0 -39 04 51.4 SCR1457-3904 12.984+0.030 3.69 • • • • • • Ml.OVI m g R 

16 13 53.5 -30 40 58.4 SCR1613-3040 12.383+0.031 3.66 • • • • • • Ml.OVI m g R 

17 56 27.9 -59 27 18.2 SCR1756-5927 12.686+0.030 3.61 • ■ • • • ■ Ml.OVI m g R 

20 18 28.7 -66 06 44.5 SCR2018-6606 12.990+0.030 — — Ml.OVI m g 

05 29 40.9 -39 50 25.6 SCR0529-3950 11.651+0.021 Ml.OVI m g 

21 01 45.6 -54 37 31.9 SCR2101-5437 12.078+0.026 3.70 ••■ ■■• Ml.OVI m- g R 

08 01 29.0 +10 43 04.2 LHS1970 13.875+0.032 3.84 9.43 374.8 esdM2.5 Ml.OVI: m- g+: Y R 8 

13 25 14.0 -21 27 06.0 LHS2734B 14.932+0.136 3.90 too far >280 Ml.OVI: m- g+: R R 

02 42 02.9 -44 30 58.7 LHS0158 9.726+0.021 3^91 6^89 191.4 M1.0[VI] m g 11 11 

17 18 35.0 -43 26 24.0 LHS0440 8.948+0.023 4.03 6.78 135.6 MIV Ml.OVI: m g+: 11 11 2 

18 22 58.7 -45 42 45.3 SCR1822-4542 12.879+0.027 3.92 ••■ ■■• Ml.OVI: m g+: R 

07 56 13.4 -67 05 20.6 WT0233 12.627+0.024 3^60 7^84 326.8 MO.OVI M2.0VI m g 6 6 6 

11 56 54.8 +26 39 56.3 LHS0318 11.797+0.018 3.65 8.77 263.8 M2.0VI: m g 18 R 

05 11 40.6 -45 01 06.0 GJ0191 5.049+0.021 3^80 7^08 161.0 sdMl.O M2.0VI m , (-0.99)"^ i Y 13 8 



Table 2 — Continued 



RA 


DEC 


Object 


Ks 


V - Ks 






Old Type 


New Type 


metallicity 


gravity 




Ref 


















km/sec 










7r 


V 


spect 




02 56 13.2 


-35 08 26.9 


LHS0162 


11.536+0.021 


3.82 


7.13 


364.1 


Ml.OVI 


M2.0VI 


m 


g+ 


6 


6 


6 




21 04 25.5 


-27 52 48.5 


LHS3620 


12.696+0.027 


3.93 


7.94 


414.6 




M2.0VI 


m 


g+ 


R 


R 






04 03 38.4 


-05 08 05.4 


LHS0186 


10.854+0.024 


4.02 


7.21 


295.6 




M2.0VI 


m— 


g 


Y 


R 






15 43 18.3 


-20 15 33.0 


LHS0406 


9.018+0.021 


4.04 


7.40 


117.4 


Ml.OV 


M2.0VI 


m— 


g 


11 


11 


2 




00 55 43.8 


-21 13 05.5 


LHS0127 


11.733+0.023 


4.06 








M2.0VI 


m— 


g+ 




R 






06 29 56.4 


-69 38 13.4 


SCR0629-6938 


12.901+0.034 










M2.0VI 


m— 


g+ 










12 24 26.8 


-04 43 36.7 


LHS0326 


11.234+0.023 


3.67 


7.69 


304.9 




M3.0VI 


m 


g 


R 


R 






15 34 27.7 


+02 16 47.5 


LHS0398 


11.502+0.025 


3.74 


7.52 


354.3 




M3.0VI 


m 


g 


Y 


19 






23 17 05.0 


-13 51 04.1 


LHS0541 


12.414+0.026 


4.05 


8.22 


420.0 


sdM3.0 


M3.0VI 


m 


g+ 


Y 


1 


8 




22 04 02.2 


-33 47 38.9 


SCR2204-3347 


11.601+0.027 


3.84 








M3.0VI 


m 


g+ 




R 






17 40 46.9 


-56 46 58.0 


SCR1740-5646 


13.195+0.040 










M3.0VI 


m 


g++ 










19 16 46.5 


—36 38 05.8 


bCR19lD— 3638 


12. 947+0. 034 


3.88 








M3.0V1 


m 


g++ 




R 






04 25 38 4 


—06 52 37 


LHS0189AB 


10 31 1+0 037 


3.94 


8.59 


128.3 




M3.0VI: 




g 


7 


7 






07 16 27.7 


+23 42 10.4 


LHS0228 


11.298+0.018 


4.20 


7.54 


300.6 


sdM2.0 


M3.0VI 


m— 


g 


Y 


R 


14 




13 42 21.2 


-35 34 50.7 


SIP1342-3534 


12.935+0.030 


4.37 








M3.0VI 


m— ; 


g 




R 






04 11 27.1 


-44 18 09.0 


WT0135 


9.834+0.020 


4.33 








M3.0VI 


m— : 


g 




5 






09 43 46.2 


-17 47 06.2 


LHS0272 


8.874+0.021 


4.29 


8.34 


86.3 


sdM3.0 


M3.0VI 


m— : 


g+ 


Y 


R 


8 




14 50 28.8 


-08 38 36.8 


LHS0381 


11.237+0.021 


3.90 


8.43 


274.9 


K7.0V 


M3.5VI 


m 


g 


Y 


19 


Y 




14 31 38.4 


-25 25 33.9 


LHS0375 


11.507+0.022 


4.12 


9.61 


157.5 


esdM4 


M3.5VI 


m 


g+ 


Y 


R 


8 




01 38 49.0 


+11 21 36.7 


LHS0144 


12.080+0.026 


4.24 


8.57 


386.6 




M3.5VI 


m— 


g 


Y 


R 






21 55 48.0 


-11 21 42.1 


LHS0515 


12.912+0.032 


4.50 


9.36 


264.5 


esdM5 


M3.5VI 


m— 


g+ 


Y 


R 


17 




05 00 15.3 


-54 06 09.0 


LEHPM3861 


13.967+0.059 


4.47 






esdM6 


M4.0VI 


m— 


g 




R 


15 




03 02 06.3 


-39 50 51.8 


LHS1490 


9.885+0.023 


4.37 








M5.0VI 


m— 


g 




R 






12 34 15.7 


+20 37 05.7 


LHS0334 


13.044+0.029 


4.98 






sdM4.5 


M6.0VI 


m 


g 




R 


17 




08 53 57.0 


-24 46 54.0 


LHS2067A 


11.571+0.023 


6.38 


9.54 


76.2 


subdwarf 


M6.0VI 


m— 


g 


R 


R 


12 





Note. — The last three columns provide references for parallaxes, V band photometry, and previous spectroscopy results. ":" indicates a questionable sub-type, metallicity, 
or gravity. All Kb magnitudes are from the 2MASS all-sky database. LHS2734A and B have zero parallax so we list them as "too far" in the table. 

''Metallicity data are from lCarnev et"ar I lll994h . 

''Metallicity data are from [Nordstrom et al. I l l2004h . 

'^Metallicity data are from lCavrel de Strobel et al. I ll200lh . There are multiple metallicities reported for Gl 12-054, so a mean value is presented. 



' ^Cavrel de Strobel et all ll200ll) . [Goldberg et alj ll2002h and lLaird et alj lll988l ) report different metallicitics, so a mean value is given. 
<=Metallicity data are from lWoolf fc WallersteinI 

iMorrison et aL I l l2003t) and I Allen et aL I ll2000l ') report different metallicities, so a mean value is given. 



References. — fl )lBessell99ol :(2')lBidelmanll985l :f3')ICavrel de Strobel et al.ll997l :( 4') ICarnev et al.ll994l:('5')ICosta fc MendeJpOOSl :f6') ICosta et al.l2005l :CplCosta et al.l2006 
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Table 3. SDSS Subdwarfs 



Object 


r 


g-r 


9-z 


Old Type 


New Type 


metallicity 


gravity 


SDSS J085843.89+511210.1 


20.00 


1.53 


3.10 


M2.0VI 


M2.0V 


m 




SDSS J093141. 85+453914.5 


19.45 


1.58 


3.08 


M2.0VI 


M2.0V 


m 




SDSS .1115117.32+011006.8 


20.35 


1.11 


3.15 


M3.0VI 


M3.0V 


111 




SDSS J003755.20-002134.2 


18.70 






Ml.OVI 


M2.0VI 


m— 


g 


SDSS J083217.77+522408.2 


19.61 


1.49 


3.14 


M2.0VI 


M2.0VI: 


m 


g+: 


SDSS J033408.64-072349.2 


20.21 


1.89 


3.35 


M2.0VI 


M3.0VI 


m 


g+: 


SDSS J090434.02+513153.9 


19.40 


1.78 


3.29 


Ml.OVI 


M3.0VI 


m 


g+: 


SDSS J161348.84+482016.0 


18.30 


1.72 


3.19 


Ml.OVI 


M3.0VI 


m 


g+: 


SDSS J112751. 35-001246.8 


20.03 






Ml.OVI 


M3.0VI 


m 


g 


SDSS J092708.10+561648.1 


19.37 


1.62 


3.02 


Ml.OVI 


M3.0VI 


m 


g+ 


SDSS J100109.54+015450.2 


19.12 


1.71 


3.25 


M2.0VI 


M3.0VI 


m 


g++ 


SDSS J101031. 13+651327.6 


19.39 


1.68 


3.30 


Ml.OVI 


M3.0VI 


m 


g++ 


SDSS J002228.00-091444.8 


18.92 


-0.89 


3.03 


Ml.OVI 


M3.0VI 


m 


g 

O 


SDSS J024501.77+003315.8 


19.38 


1.54 


3.15 


Ml.OVI 


M3.0VI 


m 


g 


SDSS J081329.95+443945.6 


19.39 


1.68 


3.35 


M2.0VI 


M3.0VI 


m 


g 

o 


SDSS J090238.75+471813.6 


19.97 


1.73 


3.21 


M2.0VI 


M3.0VI 


m 


g 

o 


SDSS J092534. 16+524442.4 


19.79 


1.71 


3.31 


M3.0VI 


M3.0VI 


m 


g 


SDSS J092745.78+582122.7 


20.47 


1.78 


3.56 


M3.0VI 


M3.0VI 


m 


g 


SDSS J095147. 77+003612.0 


18.27 


1.58 


3.15 


M2.0VI 


M3.0VI 


m 


g 


SDSS J115900.70+665214.3 


19.37 


1.61 


3.35 


M3.0VI 


M3.0VI 


m 


g 


SDSS J125919.29-025402.3 








M2.0VI 


M3.0VI 


m 


g 


SDSS J173452.52+603603.1 


18.78 


1.71 


3.36 


M2.0VI 


M3.0VI 


m 


g 


SDSS J215937.69+005536.2 


18.96 


1.61 


3.22 


M2.0VI 


M3.0VI 


m 


g 


SDSS J221500.88+005217.2 


19.08 


1.60 


3.33 


M2.0VI 


M3.0VI 


m 


g 


SDSS J145547.00+602837.3 


19.18 


1.70 


3.39 


M2.0VI 


M3.0VI 


m 


g+ 


SDSS J104320.47+010439.4 


19.16 


1.66 


3.37 


M2.0VI 


M3.0VI 


m 


g++ 


SDSS J084105.39+032109.6 


20.05 


1.45 


3.18 


M3.0VI 


M3.0VI 


m— 


g 


SDSS J091451. 98+453152.8 


19.02 


1.64 


3.49 


M3.0VI 


M3.0VI 


m— 


g 


SDSS J093024.66+554447.7 


19.14 


1.56 


2.96 


M3.0VI 


M3.0VI 


m— 


g 


SDSS J094306.37+465701.4 


19.74 


1.59 


3.29 


M3.0VI 


M3.0VI 


m— 


g 


SDSS J224854.83-091723.2 


19.85 


1.57 


3.39 


M3.0VI 


M3.0VI 


m— 


g 



- 49 - 



Table 3 — Continued 



Object 


r 


g - r 


g- z 


Old Type 


New Type 


metallicity 


gravity 


Qr>QQ T9'?c;8'?n fin mi/ii'?9 

OUOO JZOOOOU.OU — Uii^tiO.Z 


1 Q QQ 

iy .yy 




o.oo 


iVlZ.U V i 


iVlO.U V i 


ni — 


g 


ouoo J i iooui. ( ij-rUoo ( zu.d 








iVlO.U V i 


iVlO.U V i. 


ni 


g+- 


SDSS J105122. 43+603844. 8 


17.15 


1.65 


3.21 


M2.0VI 


M4.0VI 


m 


g 


QF>QQ TDQO/IOQ 7fi_LC;0'?/1 1 D 7 

ouoo juyz4zy. ( o+o/o4iu. ( 


1 8 7n 
io. ( u 


1 R7 


O.Zi 


A/ro c\\n 

iViZ . U V i 


AAA fwn 

iV14.U V i 


ni — — — 


g 








o.ov 


M9 nVT 


M4 nVT 


111 


rr 
6 


SDSS J082230.00+471645.8 


19.39 


1.67 


3.38 


M2.0VI 


M4.0VI 


m 


g 


SDSS J143930.77+033317.3 


19.38 


1.80 


3.58 


M3.0VI 


M4.0VI 


m 


g+ 


SDSS J003541. 84+003210.1 


20.15 


1.64 


3.37 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J003701. 37-003248.3 


20.15 


1.60 


3.50 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J010811. 89+003042.4 


17.34 


1.61 


3.43 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J083002.73+483251.6 


19.92 


1.90 


3.69 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J171745.22+625337.0 


18.66 


1.65 


3.57 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J221625.03-003122.5 


19.28 


1.58 


3.37 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J223802.82-082532.4 


19.71 


1.78 


3.59 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J224605.41+141640.6 


17.00 


1.57 


3.48 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J230303.49-010656.7 


18.95 


1.62 


3.42 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J235116.25-003104.8 


19.48 


1.49 


3.31 


M3.0VI 


M4.0VI 


m— 


g 


SDSS J150511. 33+620926.3 


18.61 


1.81 


3.45 


M3.0VI 


M5.0VI 


m 


g 


SDSS J230805.24+001812.7 


19.36 


1.70 


3.72 


M3.0VI 


M5.0VI 


m— 


g 



Note. — 



Column definitions are the same as in Table [2j Double lines separate each type. 
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Teff = 5500, logg^S.O 




5.0x10'^ 1.0x10^ 1.5x10'^ 2.0x10'^ 2.5x1o'* 

Wovelength (A) 



Fig. 1. — Synthetic spectra for stars having T^jf = 3500K and log g = 5.0. Red, green, 
and blue hnes represent different metaUicities, 0.0, —2.0 and —4.0. Note that the relative 
amounts of blue and red fluxes trends toward bluer objects at lower metaUicities. 
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1.0 0.8 0.6 0.4 0.2 0.0 

solor moss 



Fig. 2. — The mass-gravity relation is shown using data from Table 1 of iLopez-Morales 
(120071 ). A solid line represents a polynomial fit to all of the points, simply to be used as 
a guide. Different shades for the points represent different metallicity measurements. Stars 
with unknown metallicities are plotted as open circles. 
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0.6 



0.8 



1.0 



1.2 



Fig. 3. — The CaHl and CaH2+CaH3 versus Ti05 indices are plotted for our identified 
subdwarfs (solid circles). For comparison, known co ol dwarfs (dot s ), sub d warfs ( open trian- 
gles), and "ext r eme" subdwarfs (open boxes) from iHawley et al.l (119961 ). iGizid (119971 ). and 
Reid and Gizij (120051 ) are also shown. The dashed box indicates subdwarfs that do not 
exhibit strong spectroscopic indices, but which are either below the main-sequence on the 
HR diagram or have published metallicities [m/H] less than or equal to —0.5. A soli d line 
indicates the separation between regular subdwarfs and extreme subdwarfs adopted by lGizis 
Jl997h . 



12 



2 3 4 5 6 7 

V-Ks 



Fig. 4. — HR diagram for subdwarfs listed in Table [T] and shown in Figure [3] that have 
trigonometric parallaxes. Open circles indicate subdwarfs and filled boxes indicate "extreme" 
subdwarfs based on spectroscopic indices. Concentric circles indicate stars with CaHl index 
greater than 0.9 that are difficult to distinguish from main sequence stars at our spectral 
resolution. Open boxes indicate confirmed spectroscopic binar ies. A solid line indicates a 
fit to mai n-sequence dwarfs, primarily from iHenry et al.l (120041 ) with extra dwarf standard 
stars from I Gray et al.l (2003). The dashed line is one magnitude fainter than this solid line. 
Note that the K-type subdwarf sequence merges with the K dwarf sequence at the blue end 
of this {Mk, vsV - Ks) plot. The single point atV -Ks = 6.4 is LHS 2067A. The spectral 
types for dwarfs are given at the bottom of the figure as references. 
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6000 6200 6400 6600 6800 7000 7200 7400 
wovelength(A) 



6000 6200 6400 6600 6800 7000 7200 7400 
wovelength(A) 



Fig. 5. — Spectra for G and K-type stars from lJacoby et al.l (119841 ) (left panel, 4A resolution) 
and lSilva fc Cornelll (1992) (rig ht panel, llA reso lution) are shown. Be cause the red cutoff 



is 7400A in ljacoby et al.l ( 19841 ). the results from lSilva fc Cornelll ( 19921 ) are also plotted to 



7400A, and both sets of spectra are normalize d at 7400A. The two t ick marks indicate the 
Ba I (left) and Ha (right) absorption features. ISilva &: Cornelll (1l992l ) have removed telluric 
absorption features. 
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Fig. 6. — Top: An mid K subdwarf spectrum (black) is virtually identical to an mid K dwarf 
(red). Middle: The spectra of two early M subdwarfs differ only at CaH. Bottom: One 
subdwarf 's spectrum (black) matches an M0.5V spectral standard (red) at the blue end, but 
matches an M2.0V spectral standard (red) at the red end. The deep telluric band, O2A 
(7570A-7700A) has been removed. 
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Lil CoHl CoH2 Ti05 




6000 6500 7000 7500 8000 8500 9000 

Wavelength (A) 



Fig. 7. — Two synthetic spectra from old (red) and new (black) model grids are shown. The 
physical parameters for both models are T^jf = 3600K, log g = 5.0 and [m/H]= —2.0. Major 
differences are seen between 6500A and 7000A, and in the absorption lines of Li I (6103A) 
and Ca I (6122A and 6162A). 
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Fig. 8. — CaHl (top) and CaH2+CaH3 (bottom) indices from GAIA model grids are plotted 
against Tg//. The different line styles represent different [m/H]. The spectral types for dwarfs 
are given at the top of the figure as references. 
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6000 6500 7000 7500 8000 8500 9000 



wavelength (A) 

Fig. 9. — GAIA synthetic spectra from 4800K to 2800K are shown. All spectra are noiseless 
and have log g = 5.0. Red, green and blue lines represent [m/H] = 0.0, —1.0 and —2.0, 
respectively. Effective temperatures for each set of spectra are given above each group of 
lines. The feature marked #1 is not seen in any of our spectra. The Ca I (6162A) feature 
marked 7^2 is seen in our spectra. 
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Teff = 3500, [m/H] = -1.0 




7500 
Wovelength (A) 



Teff = 3500, logg = 5.0 




7500 
Wovelength (A) 



Fig. 10. — The top plot shows GAIA model grids at fixed metallicity ([m/H]= —1.0) and 
effective temperature (3500K). Black, red, green and yellow lines represent various log g = 
4.0, 4.5, 5.0 and 5.5, respectively. It is clear that CaH bands will be affected by changing 
gravity but Ti05 is not. The bottom plot shows model grids at fixed log g = 5.0 and T = 
3500K. Black, red and green lines represent various [m/H] = 0.0, —1.0 and —2.0, respectively. 
Note that model grids do not have telluric lines. 
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Fig. 11. — Using our spectra and metallicities measured independently by others, three 
different types of M dwarfs show the metallicity trend predicted by GA IA models. The gray 
lines represent lower metallicity stars in each pair. Metallicities from iBonfils et al.l (120051 ) 
are given in each panel for the stars. 
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Fig. 12. — K3.0[VI] to K5.0[VI] spectra are shown compared to K dwarf standard spectra 
(red). The two telluric bands (O2 A and B) have been removed for clarity. The three thick 
tick-marks represent the locations of Ca I (6162A) and Ca II (8542A and 8662A) absorption 
features. The inset plot shows the locations of the stars on the HR diagram, where filled 
circles represent subdwarfs and triangles represent main sequence standards of types K3.0V, 
K4.0V, and K5.0V (left to right). For clarity, the axis labels for the inset plot are not shown 
but are always V — Kg vs. M^s- The errors in absolute magnitudes are shown. However, 
because the errors for V — Kg are equal to or smaller than the filled circles, they are not 
shown. The solid line represents a fitted main sequence line and the dashed line is one 
magnitude fainter than the solid line. 
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Fig. 13.— Spectra of LHS 125, LHS 232 and LHS 327 (all K4.0[VI]), are shown. Note that 
the blue end of these spectra are at 6500A because of problems with these particular spectra 
between 6000A and 6500A. The two telluric O2 A and B bands have been removed. The 
triangle in the inset plot represents a K4.0V standard star. Symbols and lines have the same 
meanings as in Figure [121 
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Fig. 14.— Spectra of LHS 193A, LHS 227 and LHS 161 (all K6.0VI), are shown with our 
K5.0V and K7.0V standard spectra. The telluric O2 A band has been removed. The thick 
arrows in the HR diagram point toward sequentially higher gravity subdwarfs. Symbols and 
lines have the same meanings as in Figure [T2j 
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Fig. 15.— Spectra of LHS 164, LHS 2734A (both K7.0VI), DEN 0515-7211, and SCR0708- 
4709 (both K7.0VI:) are shown, with our K7.0V standard spectrum (represented by a triangle 
in the inset plot). The telluric O2 A band has been removed. Symbols and lines have the 
same meanings as in Figure [121 
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Fig. 16.— Spectra of LHS 165, LHS 418, LHS 244 (all MO.OVI), and LHS 300AB (MO.OVI:) 
are shown with our MO.OV standard spectrum (represented by a triangle in the inset plot). 
The telluric O2 A band has been removed. The hollow arrow indicates the shift on the HR 
diagram caused by decreasing metallicity. Symbols and lines have the same meanings as in 
Figure [H 
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Fig. 17.— Spectra of LHS 507, SCR 0709-4648, LHS 12, LHS 401 (all M0.5VI), and LHS 
521 (M0.5VI:) are shown with our M0.5V standard spectrum (represented by a triangle in 
the inset plot). The telluric O2 A band has been removed. The hollow arrow indicates the 
shift on the HR diagram caused by decreasing metallicity. The solid arrow indicates the shift 
caused by higher gravity. Symbols and lines have the same meanings as in Figure | 
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Fig. 18.— Spectra of LHS 367, LHS 299, and LHS 360, are shown (all M0.5VI:). The telluric 
O2 A band has been removed. 
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Fig. 19. — Progressions in metallicity and gravity effects in Ml. OVI subdwarf spectra are 
shown. The top set of spectra illustrates spectra with metallicities increasing from black to 

red (top to bottom, metallicity scale from m to m). The red line is our M0.5V 

standard spectrum (represented by a triangle in the inset plot). The bottom two sets of 

spectra illustrate gravity effects for metallicity scales m and m . The 

telluric O2 A band has been removed. Symbols and lines have the same meanings as in 
Figure [13 




Fig. 20. — The top set of spectra illustrates gravity effects for metallicity scale m . 

Spectra of LHS 440, LHS 1970, LHS LHS 2734B (all Ml.OVI), and LHS 158 (M1.0[VI]), are 
also shown compared individually to our Ml.OV standard. The telluric O2 A band has been 
removed. 
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Fig. 21.— Spectra for WT 233, GJ 191, LHS 406, LHS 3620, LHS 127 (all M2.0VI) and 
LHS 318 (M2.0VI:) are shown with our M2.0V standard spectrum (represented by a triangle 
in the inset plot). The telluric O2 A band has been removed. Symbols and lines have the 
same meanings as in Figure [T71 
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Fig. 22. — Spectra for six M3.0VI are shown, with our M3.0V standard spectrum in the 
top set of spectra (represented by a triangle in the inset plot). LHS 326 and LHS 228 
show metallicity changes, while LHS 326, SCR2204-3347 and SCR1916-3638 show effects 
of gravity. WT 135 and LHS 272 are M3.0VI:. Note that the K I lines at 7665A and 
7699A have appeared and blended in the O2 A band. Symbols and lines have the same 
meanings as in Figure [ITl 
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Fig. 23.— Spectra for LHS 381 and LHS 144 (both M3.5VI) are shown with our M3.5V 
standard spectrum (represented by a triangle in the inset plot). These three stars show a 
dramatic metallicity sequence for M3.5. The spectrum of LHS 228 (M3.0VI) is compared 
to that of LHS 381 to illustrate the slope difference between M3.0VI and M3.5VI (8200A- 
9000A). LHS 375 and LHS 515 have higher gravities than LHS 381 and LHS 144, respectively. 
Note that K I lines (7665A and 7699A) have appeared and blended in the O2 A band. 
Symbols and lines have the same meanings as in Figure [171 
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Fig. 24.— LEHPM 3861 (M4.0VI), LHS 1490 (M5.0VI), LHS 334 (M6.0VI), and LHS2067A 
(M6.0VI) subdwarf spectra are shown compared to our standard M3.0-7.0V standards. Only 
LHS2067A and the M6.0V standard (filled triangle) are shown in the HR diagram. Vertical 
dotted lines represent the three peaks we use to assist us to assign spectral types. The K 
I fines at 7665A and 7699A are clearly seen blended in tfie O2 A band. Symbols and lines 
fiave tfie same meanings as in Figure [TTl 
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Fig. 25. — Spectra are shown for subdwarfs from IWest et al.l (120041) s e parat ed into types 
assigned by us by matching standard spectra (red) from lBochanski et al.l (120071 ) in the region 
8200A-9000A. Other colors (black, blue, green, and yellow) represent different metallicities 
at each type. Metallicity increases from black to red. 
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Fig. 26. — These plots indicate the gravity effects for SDSS subdwarfs of type M3VI and 
M4VI. The black, blue and green lines represent increasing gravities: g, g+ and g++. 
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Fig. 27. — Observed spectra of dwarf standard stars (red lines) compared with the best 
fitting synthetic spectra (black lines). Telluric lines are not present in the synthetic spectra. 
Six regions labeled from 1 to 6 are discussed in the text. 
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Fig. 28. — Synthetic models were compared to the observed spectrum of our Ml.OV stan- 
dard. The resulting ^ values are plotted against temperature for models with five different 
metaUicities, at four different gravities each. The [m/H]= —0.5 plot has the tightest curves 
at different log g. 



- 78 - 




0,B 

0.2 I , 1 1 , I 1 , 1 1 I 1 1 1 , I 1 , 1 1 \ 1 1 , 1 I , 1 1 1 

6000 6500 7000 7500 8000 8500 9000 

wovelength (A) 




6000 6500 7000 7500 _ 8000 8500 9000 

wovelength (A) 



Fig. 29. — Comparisons of observed spectra for Ml.OV (top) and M3.0V (bottom) are shown 
against the two best fitting synthetic spectra for each. The red hne is the best fit, where 
Teff/[m/B]/log g = 3600/-0.5/4.5 for Ml.O, and 3200/-0.5/4.5 for M3.0. The blue hne is 
the second best fit, where the values are 3400/-1.0/4.0 for Ml.O and 3300/0.5/5.5 for M3.0. 
Telluric bands are marked as gray boxes. 
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Fig. 30. — The best values are shown for synthetic spectra fits to observed spectra for 
stars of types MO.OV to M5.5V (labeled relative to each type). The dashed line indicates the 
selected hmit for reliability of fitting model spectra to observed spectra — above this fine 
fits are deemed unreliable. 
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Fig. 31. — The spectrum of the subdwarf LHS 335 (black) is shown with the corresponding 
best fitting synthetic spectrum (red). 




Fig. 32. — The CaH2+CaH3 versus Ti05 indices are plotted for subdwarfs discussed in this 
paper (filled circles). For comparison, known coo l dwarfs (dots ) , sub d warfs ( open trian- 
gles), and "ext r eme" subdwarfs (open boxes) from iHawley et al.l (Il996l ). iGizid (Il997l ). and 



Reid and Gizid (120051 ) are also shown. Green circles are subdwarfs with the highest gravities 



at a given type (g with most +) in Table O The two blue circles represent LHS 400 and SCR 
1822—4542, which have higher gravities than M dwarfs. Red circles represent the lowest 
metallicity stars in Table [2] at a given type (m with most — ). Solids lines indicate pairs 
of stars having the same number of — in metallicity, but different gravities. A dotted line 
connects LHS 2734A and LHS 2734B, which is the only wide CPM binary in our sample. For 
comparison, the inset plot illustrates the indices (Ti05 vs. CaH2+CaH3) calculated from 
GAIA models for temperatures 2800-4400K. Axis labels are omitted, but both axes have 
the same ranges as the main figure. Three different metallicities, 0.0, —1.0 and —2.0, and 
two gravities, log5'=4.5 (solid line) and 5.5 (dash line), are shown. Connections are drawn 
for stars with the same temperatures but different gravities. 
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Fig. 33. — CaHl and Ti05 indices from GAIA model grids are plotted against effective 
temperature. The black, blue, red and green lines indicate log g=4:.0, 4.5, 5.0, 5.5, respec- 
tively. Note that a given CaHl and Ti05 value may correspond to many combinations of 
metallicities and gravities. 
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Fig. 34. — Spectra from iLepine et al.l (120071 ) are shown, with their pubhshed types. The 
colors (black and blue) represent different spectra, and their colors match with labels. Each 
pair has the almost the same overall spectral shape, but each star is assigned a different 
sub-type and identification as extreme or ultra. Spectra are normalized at 7500A. 
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Fig. 35. — Spectra from iLepine et al.l (120071 ) are shown for spectral types for M3.0, M5.0 
and M7.5. The black, blue and red spectra represent sdM, esdM and usdM, respectively. 
Spectra are normalized at 7500A. 



